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Abstract: It reviewed results on low-energy nuclear reaction (LENR) processes in a high-voltage (1,000-4,000 V) electric discharge 


system. The experimental results are divided into three sets: excess heat measurements; yield of impurity nuclides (nuclear ash); X-ray 


measurements. Up to 8 W of excess power was observed, with a power gain of up to 170% was seen in glow discharge experiments. Up 


to 300 W of excess power, with a power gain up to 340% was observed in experiments using a high-voltage electrolysis cell. The 


impurity nuclide yield showing a shift of up to a few per cent from natural isotopic abundances was detected by spark mass 


spectrometry, by secondary ionic mass spectrometry, and by secondary neutral mass spectrometry. X-ray emission in the range of 


0.6-6.0 keV, has been observed. Based on these experimental results we propose a phenomenological model for LENR based on the 


interaction of an electric discharge with condensed matter (of the cathode) 


Key words: Electric discharge, excess heat, impurity nuclides, X-ray emission. 


1. Introduction 


Our group has studied low energy nuclear reactions 
(LENR) in glow discharge experiments for many 
years [1]. We have identified phenomena associated 
with the Fleischmann-Pons effect in the solid cathode 
of a high-voltage glow discharge system and a 
high-voltage electrolysis cell; including excess heat 
production, transmutation, and X-ray emission. In this 
work we will review experimental results on excess 
heat production in the glow discharge system (up to 8 
W, and 170% power gain) [2], and in the electrolysis 
cell (up to 300 W, and 340% power gain) [3]. We also 
describe the results of experiments which show 
transmutation products [4]. Detailed experimental 
results on X-ray emission are discussed in another 
paper in the Ref. [5]. These include the observation of 
collimated X-ray emission from the cathode, which is a 
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new fundamental effect that has not been reported prior 
to our experiments. 

Based on these observations, we have developed 
some theoretical ideas that may be helpful in 
understanding the experimental results. We propose 
that long-lived metastable states are generated at solid 
density in the cathode with excitation energies in the 
keV range. In the glow discharge, a flux of ions with 
kinetic energies on the order of 0.5-2.0 keV is incident 
on the cathode; the appearance of X-rays in the 
0.6-6.0 keV range implies that there is a physical 
mechanism capable of converting the incident ion 
energy into long-lived (up to tens of millisecond long) 
excitation of the electron-nuclear system of the cathode. 
We observe characteristic X-ray emission of the 
cathode elements, suggesting that L and M excitation 
occurs in the cathode. Observations of excess power 
and transmutation products indicates that nuclear 
reactions occur, and it may be that these reactions occur 


as a consequence of the high effective temperature that 
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might be associated with the strongly non-equilibrium 
conditions implied by the long-lived highly excited 
states. We propose that such reactions are responsible 
for LENR effects in our experiments. 


2. Experiments 


2.1 Glow Discharge Device 


The measurements were carried out using a glow 
discharge [6] consisting of a water-cooled vacuum 
chamber, cathode and anode assemblies as shown in 
Fig. 1. The cathode holder can accept cathodes made of 
different metals, and the discharge occurs between the 
cathode and anode in the low pressure gas inside the 
chamber. Periodic pulses are supplied by a 
high-voltage power source. We used flow calorimetry, 
with independent water cooling channels supplied to 
the cathode, the anode, and to the chamber. Each 
cooling channel included two thermal sensors 
differently turned on at the input and output, and a 


cooling water flow meter. 
2.2 The High-Voltage Electrolysis Cell 


The measurements were carried out using high- 
voltage electrolysis cell; which consists of two flanges 
made of stainless still with fastening units for the anode 
and the cathode assembly as shown in Fig. 2 [7]. An 
electrolysis chamber in the form of quartz glass tube is 
set between the flanges, and another cylindrical tube 
set co-axially serves as the chamber cooling jacket. The 
electrolysis chamber volume was 95 cm’. The cathode 
units were provided with changeable cathode holders 
that allowed us to use different cathode samples with 
associated fasteners to hold the samples. The discharge 
chambers were wrapped up into a thermal insulation 
cover, witch reduced thermal losses into the 
environment. The three units of the device (the cathode, 
anode and chamber) had independent channels of water 
Each 


temperature sensors (at the input and output) and a 


cooling. cooling channel included two 


volumetric counter of the cooling water flow. The 


stirrer and the thermocouple were fixed upon the anode 
flange inside the chamber. Pd, Ni and Pt cathode 
samples were used in the experiments; the cathode 


diameter was 11 mm. 


Fig. 1 Schematic representation of glow discharge device; (1) 
discharge chamber; (2) cathode holder; (3) cathode sample; 
(4) anode; (5) cooling water; (6) X-ray output; (7) Be foil 
screens; (8) X-ray scintillator, photomultiplier detectors. (9); 
X-ray detectors different kind; (10) X-ray spectrometer; (11) 
heat insulation cover and (12) windows in heat insulation 
cover. 


Fig. 2 Schematic of a high-voltage electrolysis device, set up 
as a heat capacity and flow calorimeter; (1) electrolyte 
chamber; (2) cathode holder; (3) cathode sample; (4) 
thermal insulation cover; (5) mixer; (6) anode; (7) 
temperature sensors; (8) flow meter; (9) thermocouple and 
(10) cooling water. 
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2.3 Power Supply 


The power supply is designed to provide the glow 
discharge with periodic-pulse direct current, and 
permits the generation of desired current forms (of 
various pulse length and pulse period) to obtain the 
required current voltage operating conditions. The 
power supply consists of an autotransformer; a step-up 
transformer; a rectifier; a storage capacitor; a ballast 
resistor; and a high voltage transistor switch (Fig. 3). 
The power supply produces direct pulse-periodic 
In_ the 
experiments, we used pulse duration from 0.1 to 2.0 ms; 


current of rectangular shape of pulse. 


and a period from 4 up to 100 ms. The glow discharge 
conditions were: current (amplitude), from 30 up to 
300 mA; voltage, from 1,500 to 4,300 V. The 
electrolysis voltage in different experiments ranged 
from 500 V to 2,500 V, with the current in the range of 
0.5-2.0 A. 


2.4 Excess Heat Power Measurements in Experiments 


The excess power parameters (determined from 
input from the temperature sensors and the flow-meter), 
and the electric parameters (the glow discharge current 
and voltage), were recorded using a data acquisition 
board. The excess heat power Pey was determined as 
the difference between the output thermal power and 
input electrical power. The discharge current and voltage 


Fig. 3 Schematic of power supply; (1) autotransformer; (2) 
Step-up transformer; (3) rectifier; (4) resistor; (5) storage 
capacitor; (6) master generator; (7) guide transformer. (8) 
high-speed electronic switch based on fast IGBT transistors; 
(9) electric discharge chamber; (10) anode and (11) cathode. 


was measured with a two-channel data acquisition 
board. The use of 8-bit analog digital oscilloscope 
conversion (with 50 MHz clock frequency) allowed us 
to achieve about 1% measurement accuracy of electric 
parameters. The electric power value having been 
obtained by the following expression: 
P,, = 1/T x JU(t)(t)dt (1) 
Instantaneous current and voltage values were 
multiplied digitally and the resulting estimate for 
electric power was made; with the average electric 
power value obtained from We used a numerical 


integration based on 


1 a Pp * Dis 
P = x ) 
el T 2 


here, T (the pulse-repetition period) is set by the power 


- (rex t:) (2) 


i=0 


supply driving a generator with the accuracy no less 
than 0.1%; pi, pis representing electric power values 
for different temporal meanings (ti, ti; ). The number 
of partitions i= N was determined by the condition of 
the power correct linear approximation for each length 
ti, ti+1. The calorimetric system actual error was defined 
in tests by a tailor-made resistive heater. The 
water-cooled flow heater in size was very near to an 
electric resistance, which in its turn was close to 
discharge resistance (about 5 kQ and 1 kQ) and was 
characterized by its own small specific heat and 
adequate thermal insulation. In check experiments the 
heater was powered by the pulse-periodic power 
supply like Glow Discharge device and High Voltage 
Electrolysis Cell (HVEC). The procedure used in 
processing of results was the same as that used in 
HVEC experimentation. In so doing the electric power 
of calibration heaters was defined and measured. The 
temperature difference (in calibration) was determined 
by the formulae: 


AT, = Peta/(CwGwa) (3) 
AT¢ = Pricl(Gw Gwe) (4) 
AT ch = Peica/(CwGwen) (5) 


here, Pra, Peic, Peicn denote electric power (values) of 
calibration heaters in the anode, cathode and chamber 
water channels, respectively. 
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The heat power carried away by the anode, cathode 
and chamber cooling water channels, respectively, is 
determined in the following way: 


Pua = Cy-Gwa:AT, (6) 
Puc _ Cy -Gwc-ATc (7) 
Pucn = Cw-GwenATch (8) 


where, C,, stands for the water specific heat; Gwa, Gwe, 
Gwcn show the anode, cathode and chamber cooling 
water flow rate, respectively; AT,, ATc, ATcn denote 
the difference between the anode, cathode and chamber 
input and output cooling water temperature readings. 
The water flow rate is measured by volume flow meters 
with +0.5 % relative measurement error. 

With the availability of an excess heat power source 
(Pen) inside HVEC, its value is determined by the 
following: 

Pex = (Puc + Pua + Puch)— Pert APerror (9) 

The heat efficiency 7 (the ratio between output heat 
power and the input electrical power) was determined 
by expression 

7 = (Puc + Pua + Puch)/Pel (10) 

Testing of the system was carried out measurements 
in experiments with Pt cathode under the assumption 
that no there is excess heat power. Heat balance of the 


error information is + 3% (Figs. 2b and 12). 


3. Pre-Deuterated Pd Cathode Sample 
Preparation 


Bombardment by D2 ions in the glow discharge is 
viewed as the most effective procedure for loading 
deuterium into Pd [7]. The measurements were carried 
out using the glow discharge device consisting of a 
water-cooled vacuum chamber (1,000 cm’), and 
cathode and anode assemblies. The chamber design 
allowed the placement of cathode samples made of 
various materials on a water-cooled surface. The 
experiments were carried out using a high-current glow 
discharge in D2, with Pd samples. The power supply 
feeds the glow discharge with a periodic pulse direct 
current, and allows the generation of desired current 
forms (with various pulse length and pulse period) to 


obtain the required current voltage and current. We 
observed good Dy, loading into Pd cathodes with this 
approach. The glow discharge current was from 30 mA 
up to 300 mA, the voltage was from 500 to 1,400 V, 
and the gas pressure in the chamber was between 1 and 
25 Torr. The loading procedure can be summarized as: 

(1) First pump the gas out of the discharge chamber; 

(2) Bleed deuterium into the chamber until the 
pressure builds up to several Torr; 

(3) Switch on the glow discharge and commence 
operation. 

The Dz pressure is observed to decrease steadily 
when loading into the cathode Pd. Then the glow 
discharge is switched off, with D2 being pumped again 
into the chamber. The above procedure is repeated 
several times, the gas pressure being increased each 
time as shown in Fig. 4. During the experiment the 
D/Pd ratio was determined by summing up the pressure 
drops. 


4. Excess Power Measurements in Glow 
Discharge Experiments 


The efficiency 7 is determined in experiments under 
operational regimes in which excess heat is not 
observed. The Pd samples used in glow discharge 
calibration experiments with Xe and Kr were not 
pre-charged with deuterium, and relatively high values 
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Fig. 4 Time history of the D, pressure in the glow discharge 
chamber; (a) D, loading in Pd_ foil sample; (b) 
nanostructured Pd cathode sample, (c) Pd foil coated Au 
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for the thermal efficiency were achieved. Following 
calibration, deuterium is loaded into the Pd. In these 
experiments the current density did not exceed 
100 mA/cm’. At such values of the discharge current 
density in D2, a continuous loading of D2 into Pd ran up 
to saturation. The amount of deuterium loaded into 
palladium was determined by the volume of the gas 
absorbed from the discharge chamber. When saturation 
was achieved, the value of the D/Pd ratio was close to 1. 
Then, experiments were carried out with Pd cathode 
samples in D,, and also with deuterium pre-charged Pd 
cathode samples in Xe and Kr discharges. Relatively 
high values of excess power were achieved for 
deuterium pre-charged cathode samples in Xe and Kr 
discharges [7] (Fig. 5). 

No excess power was observed in the cathode 
samples made of pure Pd (not deuterium pre-charged) 
in Xe and Kr discharges, as in Fig. 6. The biggest 
values of power gain were observed under conditions 
when the glow discharge voltage ranged between 1,000 
and 1,300 V (Fig. 6). High values of the excess heat 
power were achieved for Pd deuterium pre-charged 
cathode samples in D2, Xe, Kr discharges. No excess 
heat power production was observed in the cathode 
samples made of pure Pd (not deuterium pre-charged) 
in Xe and Kr discharges (Fig. 7). 

We interpret this as showing two requirements for 
excess power generation: (1) deuterium should be 
loaded into the solid-state crystal lattice medium and (2) 
the crystal lattice should get an initial excitation, so that 
high-energy long-lived excited levels are created in the 
cathode sample solid. These excited conditions could 
be created by an additional source (for example by a 
flux of inert gas ions). 


5. Excess Power Measurements in the 
High-Voltage Electrolysis Cell (HVEC) 


The high-voltage electrolysis cell device was used in 
the two experimental runs: 

(1) Experiments based on heat capacity calorimeter. 
The cooling water of the cathode, anode, and discharge 


chamber, was absent in these experiments; 


(2) Experiments based on continuous flow 
calorimeter. 
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Fig. 5 Excess heat power a function of input electric power 
for a Pd cathode sample, d = 9 mm, current is 50-100 mA. (1) 
Deuterium pre-charged Pd cathode samples in D,; Xe, 
optimal (1,100-1,300 V) Glow Discharge voltage, (2) 
Deuterium pre-charged Pd cathode samples in D,, Xe. 
Non-optimal (less 1,100 V) Glow Discharge voltage; (3) 
non-deuterium pre-charged Pd cathode in Kr discharge. 
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Fig. 6 The ratio of thermal power to electrical input power 
as a function of the glow discharge voltage. Deuterium 
pre-charged Pd cathode samples in D2, Xe, Kr discharges. 
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Fig. 7 Dependence of the cathode output heat power of the 
input electric power: (1) deuterium pre-charged Pd cathode 
samples in D,, Xe and Kr discharges, d = 9 mm, optimal 
(1,100-1,300 V) glow discharge voltage and (2) non 
deuterium-charged Pd cathode, in Kr discharge. 
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Fig. 8 Dependence of the calorimeter temperature up input 
electric energy value; (1) electrolysis in H,O with the Pd 
cathode, (2) electrolysis in H,O with a Pt cathode, (3) input 
electric energy in calorimeter with Pt cathode, (4) input 
electric energy in calorimeter with Pd cathode and (5) excess 
energy with pre-loaded Pd cathode. 


5.1 Heat Capacity Calorimetry in High-Voltage 
Electrolysis Cell Experiments 


A variety of calorimetry experiments were carried 
out in the high-voltage electrolysis cell to study excess 
power generation. For control experiments, we used Pt 
cathodes (which do not load, and which produce no 
excess power). Excess power measurements were 
made on Pd cathodes that were pre-loaded with 
deuterium, and also with Ni cathodes. The electrolysis 
was done in light water. In heat capacity calorimetry, 
estimation of the excess energy is made by comparing 
the temperature increase measured with deuterated Pd 
or Ni cathodes with the temperature rise of the Pt 
control. When a Pd cathode pre-loaded with deuterium 
is used, the temperature rise is much faster.To estimate 
the excess energy, we begin with the observed 
temperature as a function of input electrical energy; for 
each observed temperature we determine how much 
input energy is needed to produce such a temperature in 
the Pt control experiment, and then determine the 
excess energy as the difference between the two lines 
(the estimated energy produced minus the input 
electrical energy) Fig. 8. The electric discharge was 
switched off when the calorimeter temperature reached 
57 °C. It is possible to develop estimates for the excess 
power 


produced from such measurements by 


determining the incremental excess energy produced 
during the time between successive temperature 
measurements. This approach was used for the 
estimated excess power in the different experiments of 
Figs. 9 and 10. We see from these measurements that 
the excess power produced reached 280 W, with a 
power gain (ratio of thermal output power to electrical 
input power) as high as 340%. 


100 -———_> 
80 Fs 


Excess heat (W) 
a 
ca) 
i 
lift 
a 


0 40 80 120 160 


200 
iC eee a 
160 
3 140 


io) 


= 120 


2 100 


(%) 


Heat power out 


0 40 80 120 160 
Time (s) 
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The big values of Excess Heat power were registered 


for different kinds Pd cathode samples (Fig. 11) and Ni 
samples (Fig. 12). 
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100 


5.2 Flow Calorimetry in High-Voltage Electrolysis Cell 
Experiments 


Our flow calorimetry in the high-voltage electrolysis 
cell is similar to the calorimetry described above for the 
glow discharge experiments. Independent cooling 
water flow was channeled to the cathode, anode and to 
the chamber. Measurements of the temperature 
difference between the input and output of the different 
channels was used to determine the thermal power 
associated with the cathode, anode and chamber 
individually. 

Experimental results are tabulated in Table 1 for Pt 
control experiments, and for experiments with Ni 
cathodes. No excess power was observed with Pt 
cathodes (Fig. 13), and the measured calorimeter 
efficiency was between 94% and 97%. Excess power 
gain was seen with Ni cathodes with good 
reproducibility; the largest excess power observed was 
230 W with a power gain of 280%. The effect of the 
mixer (stirrer) can be seen in the greater proportion of 
power showing up in the chamber flow channel with 
the mixer on. Excess heat power is produced on 


cathode (Fig. 14). 


6. Impurity Nuclides Production 


Measurements 


6.1 He Production Measurements 


Presently, the release of excess heat power is 
accounted for by “He production and the on-going 
reactions of transmutation accompanied by the 
impurity nuclides yield [8]. We had registered “He in 
the early experiments. Some of Pd cathodes together 
with the reference samples were analyzed at the 
Rockwell International Laboratory (Oliver’s group). A 
small increase in the *He concentration and a large 
increase in the “He concentration was found in the 
discharge treated Pd samples. The results of two 
experiments represented in Table 2. This is the third 
independent evidence of “He presence in the nuclear 


reaction. 


78 Experimental Results on Excess Heat Power, Impurity Nuclides and X-Ray Production in Experiments 
with a High-Voltage Electric Discharge System 


Table 1 Flow calorimetry results for high-voltage electrolysis cell. 


Pt cathode, H20 electrolyte with working mixer 


N Exp. 1 2 3 4 5 
Current* I (A) 1.13 1.6 1.1 1.4 0.64 
Voltage* U (V) 600 443 295 380 1,520 
Elec. Power P., (W) 170 178 82 132 245 
Cathode Heat power Pyc (W) 23.6 26.6 10.7 24.0 15.1 
Chamber Heat power Pucn (W) 124 128.5 57 100 145.4 
Anode Heat power Py, (W) 15.5 16.6 9.5 15.2 74.3 
Total Heat power Ps (W) 163 172 77.2 147 234 
Excess Heat power Pry (W) - - - - - 
Heat efficiency n (%) 96 96 94 95 96 


Ni cathode, H,O electrolyte without working mixer 


N Exp. 1 2 3 4 5 
Current* I (A) 0.49 0.54 0.49 0.63 0.58 
Voltage* U (V) 1,500 1,900 1,700 1,460 1,770 
Elec. Power P., (W) 183 265 207 230 212 
Cathode Heat power Puc (W) 93.2 80.3 94.2 156.4 128.8 
Chamber Heat power Pucy (W) 128.3 199 198.3 165.6 127.6 
Anode Heat power Py, (W) 57.4 81 30.6 42.7 63.5 
Total Heat power Ps (W) 279 354 323 365 320 
Excess Heat power Pgy (W) 88 89 116 137 108 
Heat efficiency n (%) 148 134 156 160 151 


Ni cathode, H,O electrolyte with working mixer 


N Exp. 1 2 3 4 5 

Current* I (A) 0.42 0.58 0.4 0.84 0.88 
Voltage* U (V) 1,110 1,070 620 654 634 
Elec. Power P.;(W) 109 151 62 138 139 
Cathode Heat power Pyc (W) 15.8 28.2 16.7 6.7 42.7 
Chamber Heat power Pucy (W) 113.5 162.5 93.8 286.5 268 
Anode Heat power Py, (W) 21 17 7.9 20.9 20.8 
Total Heat power Ps (W) 150 207.6 118 368 332 
Excess Heat power Pgy (W) 40 69 56.3 230 180 
Heat efficiency n (%) 136 137 190 280 240 


*Current and voltage values was for pulse amplitude. 


1) removal the upper 1.5 nm-thick defect | b 
6.2 Stable Impurity Nuclides Measurements (1) ; ap a Ae 
plasma etching; 


We carried out elemental and isotopic assays of 
cathode sample before and after glow discharge 
experiments to search for possible nuclear reaction 
products. This was done using X-ray fluorescent 
spectrometry, spark mass spectrometry, secondary 
ionic mass spectrometry, and secondary neutral mass 
spectrometry. In the case of secondary ion mass 
spectrometry, the procedure that we used was: 


(2) scanning the first and the second layers in 5 nm 
increments, while determining the content of the 
impurity nuclides; 

(3) removal of a layer with the thickness of 700 nm, 
and repeated scanning of the third and fourth layers in 5 
nm increments while again determining the content of 
the impurity nuclides (Fig. 15). 

The initial impurity content was a few 10 ppm. We 
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Fig. 13 The output heat power and the input electric power 
ratio for the different electric powers of HVEC discharge. Pt 
cathode sample. 
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Cathode Heat power output/EL. power input 
Fig. 14 Dependence of the output heat power and the input 


electric power ratio of the output cathode heat power and 
the HVEC input electric power ratio. (1) Ni cathode sample, 
discharge voltage is 600-1,800 V, current is 0.6-0.9 A. (2) Ni, 
500-700 V, 0.6-0.9 A. (3) nanostructure Pd, 700-1,700 V, 
0.5-0.6 A. 


Table 2. Relative content of *He and ‘He in Pd cathode 
samples from a glow discharge experiment. The cathode was 
Pd, the gas was D,, the current was 35 mA, and the 
experimental run time was 4 h. 


3 3 4 4 
# He after discharge! He initial He after discharge! He initial 
1 up to 10 time up to 100 time 
2 up to 2 time up to 35 time 


measured impurity nuclei with masses both less than, 
and greater than, that of Pd as depicted in Fig. 16. The 
main impurity elements (isotopes) with more than 1% 
content included: ’Li, a ee PN Ne, esi “Ca, Ca, 
66 Fe, Co, “Zn, Zn, BAs, Ag, Mag, og, 
"Cg, '?Cg, "Cg, |In (Table 3). The impurity 
content in the cathode sample volume was defined at 


different depths. The experiment duration in this case is 
up to 8 x 10° s. The deviation from the natural isotope 
ratio for these impurity nuclei is up to several tens of 
times. Interestingly, the following isotopes were absent: 
Ni, Ge, Ge, “Ge, Cd, and a Oi as shown in 
Fig. 17. These peculiarities are also registered within 
1 wm thick surface layer, the observed natural Pd 
isotopes ratio of the sample being changed (Fig. 18). 


7. X-Ray Measurements 


We have observed both diffuse and collimated X-ray 
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Fig. 15 Registration procedure for the impurity contents in 
the cathode samples by methods: Secondary Ion Mass 
Spectrometry (SIMS) and Secondary Neutral Mass 
Spectrometry (SNMS). (1) dirty superficial layer, (2), (3) 
analyzed layers (1), (2) scans, (4) surface of the cathode 
samples, (5) removal of a metal layer, (6), (7) analyzed layers, 
(3), (4) scans. 
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Fig. 16 Impurity nuclei content in the surface layer of the 
cathode sample (with a thickness of 100 »m) after glow 
discharge; Pd cathode, D, gas, 100 mA current, and 22 hrun 
time. 
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Table 3. Stable impurity nuclides contents on different 
cathode sample depths after experiments. 


A Impur. 1 scan 2 scan 3 scan 4 scan 
wubkde 10 nm, 50 nm, 700 nm, 800nm, 
content, % content,% content,% content, % 
hi 0.075 0.22 0.21 0.16 
Li 0.84 0.53 0.45 0.47 
NB 0.14 0.31 0.18 0.18 
a @ 0.93 0.63 0.47 0.54 
RG 0.19 0.15 0.05 0.06 
Ne 0.14 0.27 0.14 0.16 
“Ca 0.72 1.14 1.08 0.8 
“Ca 2.0 3.2 3.1 2.6 
S¢ 0.74 0.91 0.86 0.8 
OT; 0.57 0.72 0.52 0.7 
"Tj 0.25 0.14 0.31 0.14 
BTj 1.1 1.23 1.1 0.66 
Cr 0.62 0.41 0.31 0.1 
Fe 2.9 2.6 3.1 2.7 
Re 5.5 3.25 3.53 3.16 
*Co 1.0 1.0 1.4 1.5 
7H 0.21 0.43 0.54 1.0 
Ga 4.0 49 5.6 3.4 
Ge 5.1 44 5.1 6.0 
PAs 6.2 49 74 4.7 
Se 3.4 3.9 4.8 4.0 
BSe 4.5 3.45 5.8 1.4 
Br 3.0 24 2.8 - 
Se 4.0 3.4 2.5 23 
82S¢ 3.4 3.0 x 3.2 
SRb 2.2 3.4 3.3 3.6 
88S 3.1 44 4.2 6.0 
Fr 2.4 1.5 2.3 5.8 
TLE 2.8 3.0 3.0 3.4 
'2og 3.4 3.2 42 ts 
'BCq 4.0 1.8 2.8 5.1 
"4c 47 3.9 3.3 3.6 
"3th 2.2 2.5 2.3 : 


emission in the keV region in a large number of 
experiments carried out over the past decade; in this 
section we will provide a brief overview of some of our 
results [9]. Experiments were carried out using a glow 
discharge system, with a variety of X-ray diagnostics; 
including an X-ray pinhole camera, 
thermo-luminescent detectors with foils, scintillation 
detectors with a photomultiplier, and a curved mica 
spectrometer with X-ray film. Experiments were 


carried out on different metal cathode materials, 


including Al, Sc, Ti, Ni, Nb, Zr, Mo, Pd, Ta, and W. 
Different gasses were also used, including D2, H2, He, 
Kr and Xe. The power supply delivered a train of 
pulses, with pulse durations between 0.1 and 2.0 ms, 
and periods between 0.3 ms and 100 ms. We used 
currents between 30 and 300 mA, and voltages 
between 1,500 and 4,300 V; the gas pressure in the 
chamber was 3-5 torr. 


7.1 Pinhole Camera Experiments 


The X-ray emission was of sufficiently high 
intensity that it was possible to image the cathode using 
a pinhole camera. We used a transverse 0.3 T magnetic 
field to verify that the image was not due to charged 
particles; it was observed that the image was nearly the 
same with and without the magnetic field. 


7.2 Thermo-Luminescent Detector Experiments 


We used crystalline Al,O; thermo-luminescent 
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Fig. 17 Impurity nuclei ratio change in the surface layer of 
the cathode sample (with thickness of 100 pm) after glow 
discharge run; Pd cathode, D, gas, 100 mA current, and 22 h 
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Fig. 18 Pd isotope ratio change in the surface layer of the 
cathode sample (with thickness of 100 um) following a glow 
discharge run; Pd cathode, D, gas, 100 mA current, and 22 h 
run time. (1) natural ratio, (2) ratio after glow discharge. 


Experimental Results on Excess Heat Power, Impurity Nuclides and X-Ray Production in Experiments 81 
with a High-Voltage Electric Discharge System 


detectors covered with Be foils of different thicknesses 
in order to estimate the average energy of the X-ray 
emission [9]. We found that the main component was 
between 1.3 keV and 1.8 keV, with evidence of a 
higher energy component as well. These measurements 
showed that the X-ray intensity from the cathode 
surface increased exponentially with an increase in the 


discharge voltage, reaching as upper limit of 1 W. 
7.3 Scintillator and Photomultiplier Experiments 


We made use of organic with day (PMMA) 
scintillators with the optical scintillation detected using 
study the 
spatial-dependence, and _ the 


a photomultiplier and scope, to 
energy-dependence, 
time-dependence of the X-ray emission [9]. Tests with 
a transverse 0.3 T magnetic field again showed that the 
radiation was not composed of charged particles. We 
measured the transmission through 15-30 um thick Be 
foils, and the difference in transmission was consistent 
with an average energy in the range of 1.0-2.5 keV. It 
was found that the average energy for different cathode 
materials in these measurements were in good 
agreement with the thermo-luminescent detector 
above. We also did 
time-resolved measurements which showed that the 


measurements discussed 
diffuse radiation was produced while the discharge was 
on, and collimated emission was observed in bursts that 
followed the switching off of the discharge with some 
delay. 


7.4 Curved Mica Spectrometer Experiments 


Many experiments were carried out with the curved 
mica crystal spectrometer to give energy-resolved 
X-ray spectra on X-ray film [9]. A variety of spectral 
features were observed. We saw characteristic X-ray 
emission from Kr and Xe, which showed up as strong 
spectral lines on the film. We also saw weaker 
characteristic X-ray emission from transitions in the 
host metal atoms. There was a continuum apparent that 
originated from the cathodes surface, centered roughly 
near 1.5 keV, with somewhat different widths 


depending on the discharge voltage, and with an 
intensity that was weakly correlated with the cathode 
material. An example of this continuum emission is 
shown in Fig. 19. Finally, we saw extremely strong 
micro-beam emission that was sufficiently intense to 
cause a bleaching (the proper technical term is 
solarization) of the film, and which appeared on the 
film as spots (instead of as spectral lines) and curves 
(consistent with a spot which moved during the 


emission). 


8. Interpretation and Proposed Physical 
Mechanisms 


Based on the experimental results reviewed in the 
previous sections, we can consider possible physical 
processes that might be involved. Presumably, some 
excited energetic levels are formed in the cathode solid 
when its surface is exposed to bombardment by the 
ions flux generated in plasma or electrolyte medium. It 
would follow that the observed X-ray emission occurs 
as a result of de-excitation of these energetic levels. 

Deuterium ion acceleration in the glow discharge 
near-cathode region produced ion energies from 
several hundred eV up to a few thousand eV. 

Deuterium ions collide with the atoms of the crystal 
lattice ions. Presumably, some long-lived excited 
levels with energies up to several kilovolts are formed 
in the cathode solid. Under the resulting highly 


non-equilibrium conditions there exists an excited state 
21.0 
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0.6 


Specific X-ray emission intens 
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Fig. 19 Continuum X-ray emission measured with the 

curved mica crystal spectrometer from a glow discharge 

experiment with a Pd cathode and D, gas. 
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population that can be characterized by an effective 
temperature on the order of 1-3 keV, which in our view 
are the conditions necessary for low-energy nuclear 
reactions to occur. Within this framework, we can think 
about what kinds of reactions in particular might lead 
to the formation of stable transmutation products. 

The following reactions may be possible: 

(1) Pd + mD > [PdmD]* 

(2) [PdmD]*—> Pd* + *He + Heat 

(3) [PdmD]*— A* + B* > A+B + Heat 

(4) [PdmD]*— C + Heat 
[PdmD]* is a 


compound nucleus; m = 


short-lived intermediate 
1, 2, 3...; A*, B* denote 


nuclear isomers of nuclei with masses less than that of 


where 


Pd; where A, B are stable nuclides; and where C stands 
for a nuclide with a mass more than that of Pd. The first 
step in the proposed process is the formation of an 
excited state compound-nucleus. Then one of three 
possible modes is realized: 

(1) The compound nucleus may lose its excitation 
and form an exited Pd nucleus and “He; 

(2) The compound nucleus may split into two 
nuclear fragments with masses less than that of Pd. In 
so doing the two nuclei should be in excited isomer 
state (experiments show that the nuclear reactions 
energy is not produced as kinetic energy of nuclear 
fragments); 

(3) The compound nucleus may lose its excitation 
and form a stable nucleus of a heavier than Pd element. 

To determine the specific physical mechanism for 
these reactions will require additional research. One 
possible type of reaction for forming the impurity 
nuclides can be a long-range (resonant) nuclear 
reaction. The mechanism of such long-range reactions 
can be explained using as an example a specific 
transmutation reaction for Pd + D (Fig. 20a) and Pd + 
2D (Fig. 20b) [9]. The formation of significant °C 


. 93 P : : 
nuclei and “Nb nuclei was recorded in the experiments. 


In this case we assume that the reaction can proceed as 
4d + 2D > | Pd; D] > °C + Nb +7.820 MeV 
Momentum and energy conservation dictates that the 
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Fig. 20 Schematic of proposed long-ranged (resonant) 
nuclear reactions; (a) for Pd + D transmutation reaction; (b) 
for Pd + 2D transmutation reaction. 


product nucleus 13C should receive 6.8609 MeV, and 
that the product nucleus “Nb should receive 0.959 MeV 
(in connection with the reaction scheme of Fig. 20a). 
Note that a nuclear excited state (nuclear isomer) with 
an energy of 6.864 MeV (and excited level width of 6 
keV) exists for '3C, and that an excited level with an 
energy of 0.94983 MeV exists for Nb. 

The difference between the energy received by 
nuclide °C and the energy of the excited level is 3.1 
keV. If the crystalline lattice has available 1.5 keV, and 
given the 6 keV width of the excited energy level, one 
might expect that there should be a high probability for 
such a long-range (resonant) nuclear reaction to occur. 

Consider the reaction (Fig. 20b): 

Pd +2D > '[Pd2D] > Rh + SLi +11.794 MeV 

Energy and momentum conservation again dictate 
that the product 6Li nucleus should end up with an 
energy of 8.880 MeV. A nearby excited state is 
available, so a similar situation can occur as indicated 
in Fig. 20b. 

The totality of the experimental results allows us to 
assume that the energy of the excited nuclear levels of 
the product nuclei is converted into heat. The specific 
physical mechanism of such conversion will require 


additional research. 
9. Conclusions 


We have reviewed many results that we have 


obtained in glow discharge experiments and in 
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We have 
observed excess power in glow discharge experiments 


high-voltage electrolysis experiments. 
up to 10 W/cm”, and with a power gain of up to 170%; 
in the high-voltage electrolysis experiments we have 
seen excess power up to 300 W, with a power gain as 
high as 340%. These experiments clearly demonstrate 
the presence of an interesting new source of thermal 
power and energy. We have reviewed results obtained 
in our search for nuclear products, where many new 
elements and isotopes appear to have been produced. 
We have observed both collimated and diffuse X-ray 
emission in the keV regime, which constitutes a new 
fundamental effect that we discovered. Finally, we 
have presented a discussion of our results in terms of 
conjectures and hypotheses concerning physical 
mechanisms and reactions that may be involved (but 
which will require further study in order to determine 


whether they are correct or not). 
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Abstract: Since 1990s, the use of deicing salts (i.e., chlorides) has dramatically increased in areas with heavy snowfall in Japan. As a 


result, the water mixed with salts has accelerated the damage of the reinforced concrete (hereafter, RC) structures. Recently conducted 


inspection results of RC bridges have reported that many of the damages or deteriorations are observed at the girder ends and abutments. 


This is caused from the water leakage due to the aged expansion joints. In general, the cost for repairing the damaged RC structures is 


much higher than that for renewing the expansion joints. Therefore, to prevent these damages, we developed a new highly durable 


jointless system, named RC plug joint, for existing RC bridges with a bridge length less than 40 meters. The RC plug joint connects the 


abutment’s backwall to the superstructure’s deck using reinforcing steel bars and fiber reinforced concrete. The newly developed RC 


plug joint system can prevent water leakage and allow for a smooth ride of vehicles at the joint. This paper will explore and discuss the 


development of the RC plug joint, analysis of conducted investigations, and future installation methods. 


Key words: Jointless system, RC plug joint, field investigation, leakage assessment, travelling performance. 


1. Introduction 


The use of deicing salts dramatically increased in 
heavy snowfall areas in Japan as a result of studded 
tired being banned in 1990 [1]. Fig. 1 shows the 
volume of deicing salts used on the Hokuriku 
Expressway. The volume of deicing salts used has 
steadily increased. The aged bridge expansion joints 
with the deteriorated waterproof function allow the 
remnants of dissolved deicing salts, rainwater and 
melted snow to drop down to the girder. As a result, 
many reinforced concrete (hereafter, RC) girders, 
abutments and piers have been damaged by salt attack. 
Figs. 2 and 3 show the typical damage of an RC girders 
end and abutments due to leakage from aged expansion 
joints. The delamination of concrete and exposure of 
reinforcing bars (hereafter, rebars) are found in the locally 
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deteriorated concrete girders and abutments. Therefore, 
we developed a new highly durable jointless system, 
named RC plug joint, for existing RC bridges with a 
bridge length less than 40 m. The RC plug joint connects 
the abutment’s backwall to the superstructure’s deck 
using reinforcing steel bars and fiber reinforced concrete. 
This paper will explore and discuss the development 
structural detail, 
conducted investigations. 


process, installation method and 


2. Investigations of the Aged Expansion 
Joints 


2.1 Investigation of Water Leakage 


Fig. 4 shows the results of the investigation of water 
leakage from several types of the expansion joints in 
the Hokuriku Expressway. Water leakage was reported 
in all kinds of existing expansion joints. It was also 
reported that salt attack occurred in RC structures 
around the aged joints [2]. The sealed comb joint was 
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Fig. 1 The volume of deicing salts in the Hokuriku 
Expressway. 


Fig. 3 Typical damage of an bridge abutment. 
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Fig. 4 Investigation of water leakage from several types of 


expansion joints [2]. 
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expected to be the least leaking type among the 
conventional expansion joints. However, its leakage 
reached around 30%. Much of the leakage was caused 
by accumulated dust and sand within the drainage. The 
accumulated dust and sand in the drainage created a 
narrow path allowing rainwater to overflow. In 
addition, the deicing salts corroded the drainage 
channel, thus increasing the leakage. Fig. 5 shows 
water leakage from aged expansion joints. 

In the past, we worked to replace or improve the 
aged expansion joints, but it proved too difficult to 


prevent the water leakage completely. 
2.2 Typical Damage of Expansion Joints 


As a result of heavy traffic, the bridge expansion 
joints were found to be the most damaged part in the 
bridge elements. In addition, the heavy cyclic traffic 
loads caused poor conditions around the expansion 
joints. 

Fig. 6 shows a typical damaged area around the 
aged expansion joints. The remnants of dissolved 
deicing salts, rainwater and melted snow eventually 
saturate the joints and leakage begins. 

Fig. 7 shows the accumulated amount of the local 
damages and repairs due to the leakage from the aged 
expansion joints. The locally deteriorated areas were 
inspected between 2 m from the edge of RC bridge 
girders and abutments. The evaluation of the inspection 
data shows the number of the locally deteriorated areas 


Fig. 5 Water leakage from the sealed comb joint. This is 
one of the most serious problems on managing expressway 
bridges. 
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Fig.6 Leakage from an aged expansion joint and damaged 
concrete area was occurred. 
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Fig. 7 Results of the inspection of bridges in the Hokuriku 
Expressway. 


has dramatically increased since 2000. Currently, the 
percentage of the repairs (i.e., total repaired area 
divided by total deteriorated area) is around 40%, 
forcing the continuation of the repair work inevitable. 
In general, the repairing cost of structural damages by 
salt attack is much higher than the cost for renewing the 


expansion joints [3]. 
3. Outline of Reinforced Concrete Plug Joint 
3.1 Basic Concepts 


Fig. 8 shows the side view of the RC plug joint. In 
this system, the aged expansion joint is removed 
completely, and fiber reinforced concrete fills the 
space. Then, a waterproof membrane is applied on the 
surface of the fiber reinforced concrete. Lastly, a new 
asphalt pavement layer is completed. 


3.2 Structural Characteristics 


The RC plug joint connects the abutment’s backwall 


(i.e., parapet) to the superstructure’s deck using fiber 
reinforced concrete. This jointless system is applicable 
to existing concrete bridges whose length is less than 
40 meters. Due to temperature change, deformation 
occurs as shown in Fig. 9. However, it was confirmed 
that the strain in the rebars was low level. As a result, 
the prevention of leakage in addition to creating an 
easier passage for vehicles at the expansion joint is 
achieved. Therefore, the superiority of this system 
compared to the previously used expansion joint 
devices 1s quite clear. 


3.3 Details of RC Plug Joint 


A detailed view of an RC plug joint is shown in 
Fig. 10. The RC plug joint completes the gap between 
the girder and abutment with fiber reinforced concrete 
for the main purpose of preventing water leakage from 
the expansion joint. In addition, the deck waterproofing 
and the road asphalt pavement are applied on top of the 
RC plug joint. This double protection system 
guarantees increased waterproof protection in 


comparison with the previously used type of expansion 
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Fig. 8 Side view of the RC plug joint. 
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Fig. 9 Deformation of the abutment and foundation. 
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Fig. 10 Details of the RC plug joint (Unit: mm). 


joint system. A concrete adhesive material is applied to 
the interface of the concrete to ensure bonding and 
cohesion of both the old and new concrete with the 
steel anchor bars. Furthermore, it is predicted that the 
RC plug joint could possibly collapse in the event of a 
large-scale because the 


earthquake length of 


embedding rebar is 6 times the rebar diameter 9. 
3.4 Installation of RC Plug Joint 


The installation of the RC plug joint is executed in 
phases, to avoid traffic interruption. The repairs are 
usually conducted in daytime. 

Table | shows the mix proportion of fiber reinforced 
concrete. The repairing concrete has short cure time 
and compressive strength reach to 24 N/mm within 
only 3 h. Moreover, the tensile strength of fiber 
reinforced concrete is stronger than conventional 
concrete. 

In many cases, the replacement of old expansion 
joints in expressway bridges has weakened the existing 
concrete surface and corroded the rebars of these 
structures. In addition, the durability of RC plug joints 
is greatly affected by the adhesion strength to the 


Table 1 Mix proportions of fiber reinforced concrete. 


Filled material 


Superstructure ® :The rebar diameter 


existing concrete. For this reason, we carefully 
removed the existing aged joints, the concrete 
delaminated areas and the remaining loose materials. 
We also conducted a visual inspection and hammer 
tapping test to assess the soundness of the existing 
concrete surface. Then we coated the concrete surface 
with the adhesion material and poured concrete. Fig. 11 
shows the installation of the RC plug joint. All 
processes were completed within 8 h. During the 
following days, we set a waterproofing membrane on 
the concrete surface and paved the asphalt layer. 


4. Results of Field Investigation 


4.1 Leakage Assessment 


As of 2012, we have replaced the existing joints with 
the RC plug joints at 85 locations within the area under 
supervision of Kanazawa Branch, Central Nippon 
Expressway Co., Ltd. Fig. 12 shows the number of 
installed RC plug joints and results of leakage 
assessment. Unfortunately, the water leakage from the 
RC plug joints was found at three locations. The faults 
of RC plug joints occurred due to the debonding of the 


; ; Weight per volume (kg/m*) Admixture 
Fib terial WIC (% /a (Yo 

iber materia (%) s/a(%) W C S G Fiber (C x%) 
Polypropylene fiber 39.7 51.9 174 438 864 807 11.6 2.0 
Vinylon fiber 40.3 52.6 184 457 853 775 16.5 2.0 
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Fig. 12 Leakage assessment of RC plug joints. 


existing surface concrete. After much discussion over 
this fault, we decided to use an adhesion material (i.e., 
epoxy resin) on the existing surface concrete for future 
repairs. 

There has been no subsequent pavement damage 
such as potholes, cracks or leakage from the RC plug 
joints. From these findings, it is confirmed that the 
reliability of the RC plug joint proves thus far 
successful and we intend to continue monitoring their 


performance. 
4.2 Travelling Performance 


In order to assess the effect of vibration reduction 


(e) (f) 
Fig. 11 Installation of the RC plug joint: (a) removing the aged joint; (b) setting rebars; (c) painting with the adhesion 


material; (d) pouring concrete; (e) curing the concrete and (f) paving the asphalt layer. 


due to traffic loading after the installation of an RC 
plug joint, the vibration response of both, the bridge 
and test car was measured. Fig. 13 shows the 
reduction effect of the bridge response before and 
after the installation of the RC plug joint. These 
results show that the RC plug joint has a great 
reducing effect on the bridge’s vibration. After 
comparing the results before and after the RC plug 
joint was installed, it is observed that a smoother 


passage for vehicles was also achieved. 
4.3 Measurement of Strain of Anchor Bars 


Fig. 14 shows the measured strain of anchor bars 
under traffic loads. Fig. 15 shows the strain of anchor 
bars due to seasonal temperature variation. It was 
concluded that the strain due to temperature variation 
was 50 times higher than that caused by the traffic 
loading. As a result of this discrepancy, future attention 
must be paid to the temperature variation during the 
design phase. Furthermore, a structural analysis on the 
effect of temperature variation using a frame model 
concluded that the permissible expansion and 
contraction of the RC plug joint should be limited to 
20 mm for bridge lengths less than 40 m. 
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Fig. 13 Velocity response of the bridge caused by test car. 
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Fig. 14 Measured strain induced by test cars. 


5. Conclusions 


This paper discusses the development of the RC 
plug joint system as well as analyzes the results of the 
investigations conducted. This particular method of 
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Fig. 15 Measured strain due to temperature variation. 


installation is expected to be the most cost-effective 
one. In the next phase of development, we propose to 
work on a simple design and establish a standard of 
installation, while conducting further replacement of 
existing deteriorated expansion joints. Although the 
safety, serviceability and durability of the RC plug 
joints have been confirmed, the authors intend to 
conduct further field surveys and monitor the 


long-term performance of the RC plug joints. 
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Abstract: The possibility of developing a stellarator-based neutron source designed for the nuclear reaction initiation in the blanket of 


hybrid reactor is studied. An analog of the Large Helical Device (LHD) stellarator design, with linear dimensions increased by a factor 


of 1.5 is taken for the magnetic system. Plasma parameters and the deuterium-tritium (DT) mixture fusion power are calculated using 


the space-time numerical code under the assumption of the neoclassical transport in the ambipolarity regime. Using the 10 MW plasma 


heating sources, it is possible to obtain the DT fusion power of one-to-two tens MW. 


Key words: Stellarator, subcritical reactor, ambipolar electric field, neoclassical transport, neutron source, fusion power. 


1. Introduction 


Investigations into the controlled thermonuclear 
fusion are going on during more than fifty years, but 
there is no last answer to the problem. Putting into 
operation of an experimental fusion reactor ITER, 
planned by the end of the current decade, will be only 
the next, though important, step on the way of problem 
solving. 

At the same time, there are existing long-discussed 
alternative paths for development of power reactors 
using the nuclear fusion as a neutron source for power 
multiplying by the nuclear fission of heavy elements or 
for the accumulation of a fuel to be used in nuclear 
reactors. In the first case for energy multiplying one 
can use uranium enriched by isotope **°U having the 
multiplication factor less than unity. In another case the 
fusion-produced neutrons participate in nuclear 
transformations of nonfissionable heavy elements into 
elements which can be used as a fuel for nuclear power 


plants. For this purpose, it is possible to use 
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uranium-238 and thorium-232. These elements absorb 
neutrons and after subsequent decays are transformed 
into plutonium-239 and uranium-233 respectively. 

In one DT-fusion event the energy of 17.5 MeV is 
released with a-particle and neutron production, and 
one thermonuclear neutron fission of a heavy element 
gives near 200 MeV. And neutrons, produced as a 
result of nuclear fission, can evoke new fissions. A 
number of such fissions is determined by the 
multiplication factor K,,, which at K,, < 1, according to 
the Ref. [1], equals to 7 = 1/(1 — K,,)/v% where, visa 
number of neutrons produced by one fission event. In 
the case of K,, > 1 the system comes into the 
uncontrolled chain reaction regime. Thus, when, for 
example, the multiplication system provides the 
number of fission events per one external neutron 7 = 
10, the reactor power can be increased more than by a 
factor of 100 as compared with the power of a 
thermonuclear neutron source. Consequently, to 
develop a Giga-Watt-power reactor, it is sufficient to 
have a thermonuclear device with a fusion power of 
10 MW. 

Unlike a traditional reactor based on the nuclear 
fission, where the multiplication factor is equal to unit, 
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in the reaction under consideration the probability of 
uncontrolled chain reaction appearance at K,, < 1 takes 
no place. 

Combination of fusion reactions and nuclear fission 
in one and the same device has determined the name of 
such systems, as hybrid reactors. Fusion-fission hybrid 
reactors have been considered as power sources as 
early as investigations into the problem of controlled 
thermonuclear fusion commenced. A detailed analysis 
of such reactors, developed on the base of 
thermonuclear traps -tokamaks, has been carried out in 
the USSR by Golovin and Orlov [2, 3]. A design of a 
Hybrid Thermonuclear Tokamak Reactor (HTTR) has 
been proposed Velikhov [4]. On the initiative by IN. 
Golovin the Kharkov researchers started designing of 
an experimental fusion facility TB-0 [5]. Unfortunately, 
the project was not realized. 

An accident at the Chernobyl nuclear power plant 
has reduced for a certain time the interest in hybrid 
systems and livened up the works in the direction of 
developing clean fusion reactors, in particular, the 
International Tokamak Exsperimental Reactor (ITER) 
project. However, to solve the energy problems by 
means of clean fusion reactors is a long-term outlook, 
therefore, during the recent years the interest in 
developing hybrid systems revived [6-9]. Feasibility of 
these designs become real as at modern thermonuclear 
facilities a nuclear fusion power of 10 MW has been 
reached [10, 11]. 

Up to date fusion-fission hybrid reactor designs were 
developed, in general, on the base of tokamaks. A 
tokamak’s disadvantage is that to maintain the 
stationary fusion reaction in the system is impossible. 
Therefore, there are some technical difficulties for the 
development of a stationary power source on the base 
of tokamaks. There is no such a problem in a stellarator, 
where to maintain the fusion reaction is possible. At the 
same time, experimental studies on the fusion plasma 
in stellarators have not yet reached the level of research 
though — the 


maintenance observed in experiments obey nearly the 


at tokamaks, features of plasma 


same laws as these noticed in tokamaks [12]. 

Theoretical investigations into the plasma loss 
mechanisms in the stellarator systems resulted in the 
development of a neoclassical transport theory [13-15], 
permitting to take into account an ambipolar electric 
filed leading to the flattening between the ion and 
electron plasma fluxes. In Refs. [16-18], a comparison 
is done for experimentally measured electric field 
profiles in the Wendenstein VH-AS stellarator (W7-AS) 
and these calculated by the neoclassical theory. In the 
calculations of E, the plasma density and temperature 
profiles measured experimentally were used. A 
satisfactory accordance between the theory and 
experiment is obtained. 

In the paper by Wakasa et al. [19], experimental 
results obtained in three stellarators (LHD, W7-AS and 
Spain stellarator TJ-II) are compared with neoclassical 
theory data. The electric field values were calculated 
using the Drift Kinetic Equation Solver (DKES) code. 
A good agreement between the calculated and 
experimental results is shown for LHD and T)J-II 
facilities, at least for the major plasma part. A certain 
correspondence between the experiment and 
neoclassical theory for the central region of the plasma 
column was obtained for the LHD and TJ-II facilities 
by Olaya [20]. 

Results of the comparison between the plasma 
confinement theory and experiments in stellarators 
give grounds for forecasting the expected plasma 
parameters in experimental facilities of the next 
generation, in particular, in the systems that might be 
used for the development of neutron sources. The 
transport coefficients in the neoclassical theory are the 
functions of local plasma parameters that enable to 
obtain their spatial distributions by solving the set of 
equations such as diffusion equations and heat 
conduction equations. Similar solutions were obtained 
in the calculations of plasma parameter for different 
variants of the stellarator-reactor: experimental [21], 
demonstration [22] and power [23] reactors. It has been 


shown that under certain assumptions on the processes 
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of fuel 
stellarator-reactor, with comparable dimensions, can 
these of 


loading and plasma _ heating, the 


have energy characteristics like to 
corresponding-type tokamak-reactors. 

The present work is aimed to the study of a 
possibility of developing the stellarator -based neutron 
source designed for the initiation of nuclear reactions in 
the hybrid reactor blanket. An analog of the LHD 
stellarator design was taken as a magnetic system with 
linear dimensions increased by a factor of 1.5. Thus the 
stellarator parameters become comparable with the 
parameters of existing experimental facilities of a 
tokamak type. In addition, the sizes of the system 
provide enough space to accommodate the breeding 
blanket around the plasma. 

The plasma parameters and the DT mixture fusion 
power were calculated by the space-time numerical 
code earlier used for the reactor parameter calculations 
in Refs. [21-23]. 


2. Set of Equation and Numerical Model 


The set of equations solved in Refs. [22, 23] 
describes the space-time behavior of the plasma in the 


stellarator reactor. The system assumes a 


one-dimensional spatial dependence of plasma 


parameters on the averaged minor plasma radius. 
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The first two equations describe the space-time 
behavior of the electron Eq. (1) and the ion Eq. (2) 
temperature. The third one is a plasma particle 
diffusion equation. The first term on the right side of 
equation for the electron heat conduction Eq. (1) 


determines the temperature variation as a result of the 
heat flux corresponding to the neoclassical theory 
Kovrizhnykh [14]. The second term on the right 
describes the heating of electrons produced in the 
process of fusion by 3.52 MeV a-particles. An equal 
content of deuterium and tritium in the reactor plasma 
is assumed. The coefficient Ky determines the part of 
a-particle energy given to electrons. It is supposed 
that 95% of a-particle energy is absorbed by the 
plasma. The rest terms on the right sides of Eqs. (1) 
and Qn. are the 


heating of ions and electrons from the external sources, 


and (2) denote respectively: Qh; 


Q.-electron-ion heat exchange as a result of Coulomb 
collisions, Q,-bremsstrahlung, Q,-cyclotron radiation, 
QOr-energy change of particles in the ambipolar electric 
field during their radial motion. The cyclotron radiation 
is well-reflected from the chamber wall, therefore in 
the calculations the reflection coefficient for this 
radiation is taken equal to 0.95. Consequently, only 5% 
of the power is absorbed by the reactor walls. 
Calculations were carried out with the use of 
different models of plasma heating by the external 
sources: separated heating of electrons and ions or 
simultaneous heating of both components. Also, it has 
been assumed that the specific heating power is 
proportional to the plasma density: Qr; = gajN, where j 
denotes electrons or ions according to the equation. 
The total heating power was determined by the integral 


over the plasma volume: 


Phy = $, Qnjdv (4) 


The specific energy release due to the effect of the 
ambipolar electric field has been taken by the value 
QO; = -S;E, under conditions when the diffusion fluxes 
of electron and ions are equal: S, = S;. 

Here, as in the Ref. [23], the heat fluxes and 
diffusion particle fluxes determined in the Ref. [14] are 
used. It has been assumed that the plasma electrons are 
in the regime when the diffusion coefficient has an 
inverse relationship with the collision frequency v, and 
the ion diffusion obeys to the relationship Vv . 
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The density equation Eq. (3) contains on the right 
side, besides a diffusion term, a term with the particle 
source Ss , used for the simulation of the plasma 
density maintenance at an approximately constant 
level. Diffusion plasma losses were compensated by 
the model of the fuel pellet injection. It has been shown 
earlier [23] that the fuel injection into the central 
plasma column region is optimum one from the point 
of view of plasma confinement. 

The set of Eqs. (1)-(3) was enlarged by boundary and 
initial conditions. Besides, the spatial derivatives of 
density, temperature and potential were assumed to be 


zero in the plasma column center: 
eee, =0, N.=0 =0, Os =0, (5) 
and the initial conditions satisfy the expressions: 


Tyo =Ty0l-¥) +5, Nyao=NI-P)+5, (6) 


jt=0 


where j denotes the particle sort, x is the current radius 
normalized to unit. The constants 6, and dy determine 
the values of the plasma density and temperature at the 
plasma column edge. The index /, being an integer, in 
the most of calculations was assumed as / = 2. 

The ambipolar electric field value was determined 
from the equation for the equality between the ion and 
electron fluxes S, = S; at each step of the numerical 
code space network. In the present paper the electric 
field in every case was determined by the left (ion) root 
of the equation for the equality of fluxes giving its 
negative values. Such a choice is not optimum from the 
point of view of reactor confinement properties. A 
regime in which positive values of the electric field 
(electron root) are realized [14] might be more suitable. 
The electron root realization occurs, as a rule, with 
lower values of fluxes S;. However, the solutions with 
the right root do not exist in every point of the plasma 


radius, while the negative root is always presented. 
3. Calculation Results 


In the neoclassical theory the heat and diffusion 


fluxes of plasma electrons depend on the effective 


helical ripple €, . Its value in different experimental 
facilities is varying within the wide ranges, and even 
under conditions of a single facility the value of €, 
undergoes appreciable changes when the vertical 
magnetic field and the plasma pressure / are 
varying [24]. In this paper, as in the Ref. [23] for 
simulation of the amplitude of magnetic field helical 


disturbances the following expression is used: 
0,87 
Ep = Ep +0,05r (7) 


where, & o = 0.01, that gives at the plasma edge 
Ehmax = 9.06n. In the LHD stellarator with magnetic 
axis shifted by the radius R = 3.53 m, the effective 
ripple at the plasma edge is equal to 0.05 [25]. 

Plasma density maintenance at an approximately 
constant level was simulated by the fuel pellet injection 
under the assumption of the triangular model ablation 
on = Np (1 


injection center and A is the half-width of injection 


| a2 y| /A), where a, is the radius of 


region. The pellet size was equal to 1% of the total 
number of particles in the reactor plasma. When the 
complete particle number decreases below 0.99 the 
next pellet is injected. The energy loss by heating the 
injected particles is taken into account via the condition 
of thermal plasma energy conservation. The energy 
loss by pellet evaporation was not taken into account 
because of its small quantity. 

As for geometric parameters, an analog of the LHD 
stellarator design, with linear dimensions increased by 
a factor of 1.5, was taken with the minor plasma radius 
of one meter (a = 1 m) and the major tore radius 
R= 6 m. In all the calculations a basic magnetic field 
was taken equal to 5 T. 

The problem was reduced to the calculation of the 
space-time behavior of plasma parameters with the 
given levels of the plasma heating power and in the 
regime of fuel pellet injection. In the most of 
calculations a 10 MW power from the external source 
was loaded into the plasma. In a different of 
calculations the power was loaded wholly into ions or 
separately into electrons. Also, there were calculations 


94 The Parameters of the Stellarator as a Neutron Source for a Subcritical Reactor 


on the simultaneous heating of ions and electrons with 
equal parts of contribution. 

It has been shown in Ref. [21] that the maximum 
plasma parameters, from the point of view of heating 
and confinement, are reached in the case of fuel 
injection into the plasma column center. Therefore, the 
most of calculations in the present study were 
performed in the variant of fuel pellet injection into the 
plasma center. According to the calculation results the 
fusion power levels of about 10 MW are reached in the 
system under consideration at an average plasma 
density value of 7 x (0-8 x 10" mi, Fig. 1 presents 
the time dependences of the DT-fusion power, when 
the power from the external heating sources was loaded 
by equal portions into ions and electrons, for two 
average plasma density values. In this case the 
injection has been carried out into the plasma center 
with a half width equal to the plasma column radius (A 
= 1). The heating process duration in this calculation 
was 28.2 s. The maximum fusion power for the plasma 
density of 7 x 10’? m® was equal to 11 MW and 13.5 
MW in the second case. 

The temporal behavior of average values of the 
electron and ion temperature, under conditions given in 
Fig. 1, is shown in Fig. 2. The electron temperature is 
slightly higher than the ion temperature. It can be 
explained by the influence of two processes: the 
a-particle energy absorption by electrons and the 
energy exchange between ions and electrons when the 
in the 


ambipolar electric field having a negative value 


potential difference has been overcome 


practically throughout the plasma column cross-section 
(Fig. 3). Note, that at the lower value of the plasma 
density its temperature value is higher. 

Despite the fact that the pellet evaporation region has 
a width equal to the plasma radius, the plasma density 
distribution is very peaked. A maximum value of the 
density in the plasma column center reaches 107 m®, 
while in the peripheral part of the plasma radius the 
density value is significantly decreased. 

As the average value of the plasma density increases 


its maximum values in the plasma center become 
higher. At the same time the profiles of the ambipolar 
electric field value practically have no difference. Note, 
also, if the average plasma temperature values are less 
than 4 keV the negative values of the ambipolar electric 
field exceed 20 kV/m. 

Fig. 4 illustrates the temporal behavior of the 
bremsstrahlung power and cyclotron radiation power 
for the case of the plasma density <n> = 7-10" m°. 
Note that in the data of Figs. 1-4, the impurity influence 
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Fig. 1 Total DT-fusion power as a function of the time for 

two values of the plasma density, 7-10°° and 8-10°° m®°, in the 

case of fuel injection into the plasma center with evaporation 

half-width A = 1 and source heating power P); = Phe = 5 MW. 
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Fig. 3 Radial distributions of the plasma density and 
electric field for their two average values at the end of 
heating process (P); = Phe = 5 MW, A = 1, t = 28.2 s). 


has not been taken into account. The total power of 
losses, caused by the bremsstrahlung and cyclotron 
radiation, amounts to 6 MW that is about 50% of the 
total power released in the plasma (10 MW is the 
heating from external sources and about 2 MW from 
a-particles). 

Fig. 5 shows the radial distributions of plasma ion 
and electron temperatures for two values of the 
plasma density: 7 x 10’ and 8 x 10’. In the central 
region of the plasma column the temperatures of ions 
and electrons practically are equal. In the peripheral 
part of the plasma the difference in the temperature 
values becomes sufficiently great. Most likely in this 
region the collisional heat exchange between 
components is insignificant because of a low density. 
The difference in the temperature values determines 
the energy exchange as a result of overcoming the 
potential difference created by the ambipolar electric 
field. The fusion powers and the plasma component 
temperatures as a function of the heating time at 
different sizes of the pellet evaporation region A are 
shown in Figs. 6 and 7. 

The nuclear fusion power in the case of a small 
ablation region width (A = 0.5) at a certain instant of 
60 MW. 
realignment leads to the fusion power decrease to about 


time reaches Subsequent parameter 


40 MW. Corresponding changes in the plasma 


temperature are observed in the time-temperature 
curves (Fig. 7). The ablation width increase leads to the 
plasma temperature and fusion power decrease. The 
difference between the temperatures of ions and 
electrons in each of the above-mentioned cases does 
not exceed 100 eV. 

The results given in Figs.1-7 were obtained by the 
fuel pellet injection into the plasma column center with 
different values of the evaporation region width. 
The calculations show that the fuel injection into the 
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Fig. 4 Bremsstrahlung and cyclotron radiation power as a 
function of the heating time for the plasma density value 
<n> = 7 x 10’° m® in the case of fuel injection into the 


plasma center with evaporation half-width A = 1 and source 
heating power P); = Phe =5 MW, Z=1. 
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Fig. 5 Radial distributions of the plasma temperature for 
two average plasma density values at the end of the heating 
process (P),; = Phe = 5 MW, A = 1, t = 28.2 s., Z = 1). 
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Fig. 6 DT-fusion power in the plasma heating process as a 
function of the time with different half-widths of the fuel 


pellet evaporation region A = 1, 0.75, 0.5. <n> = 0.7 x 10” m®, 


Phi = Phe = 5 MW, Z = 1. 
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Fig. 7 Plasma component temperature as a function of the 
time in the plasma heating process with different 
half-widths of the fuel pellet evaporation region A = 1, 0.75, 
0.5. <n> = 0.7 x 107° m®, P,,; =P,. = 5 MW, Z = 1. 


peripheral plasma region leads to the significant 
change in the spatial distribution of plasma parameters. 

Figs. 8-10 demonstrate a case of fuel injection into 
the region of the plasma column half-radius with a 
half-width of the pellet evaporation region A = 0.24. In 
this case a maximum of the plasma density is formed in 
the fuel injection region with a deepening in the plasma 
column center. Maximum negative values of the 
ambipolar electric field are shifted relatively to the 


density maximum towards the plasma edge. 
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Fig. 8 Radial distributions of the plasma density and 
ambipolar electric field in the case of fuel injection into the 
region of the plasma column radius half-width (x; = 0.5, P); = 
Pie = 5 MW, A = 0.24, <n> = 7 x 10" m®, t = 36.2 s, Z = 1). 


The radial distributions of ion and _ electron 
temperatures are shown in Fig. 9. In the plasma center 
the temperatures have a maximum, unlike the plasma 
density profile the maximum of which takes place in 
the fuel injection region. Here the ion temperature is 
slightly higher than the electron temperature. The 
temperature equalization takes place in the fuel 
injection region. Near the plasma edge the electron 
temperature is higher than the ion temperature. It is due 
to the energy exchange between the components in the 
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Fig.9 Radial distribution of the plasma temperature in the 
case of fuel injection into the region of the plasma column 
radius half-width (x; = 0.5, Py; = Phe = 5 MW, A = 0.24, <n> = 
7x 10" m°, t= 36.2 s, Z = 1). 
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presence of the radial electric field. 

Fig. 10 shows the distributions of bremsstrahlung 
and cyclotron radiation specific powers. Maxima of 
radiations correspond to the plasma density maximum. 
In the regions with a low plasma density the cyclotron 
radiation power is higher than the bremsstrahlung 
power. The ratio between the radiation power values 
changes to the opposite one in the region with a higher 
plasma density. Note that the nuclear fusion 
power in this variant of fuel injection is at a level of 
0.1 MW. 

The above-given calculation results were obtained 
under the assumption of the impurity absence in the 
plasma (Z = 1). Nevertheless, the bremsstrahlung 
withdraws a significant part (to 30%) of the fusion 
power and external heating sources released in the 
(Figs. 1-4). At the 
bremsstrahlung power strongly depends on the value of 


plasma same time, the 
the plasma ion charge Z. Fig. 11 illustrates the temporal 
behavior of the fusion power for the cases Z = 1, 1.2 
and 1.5. The calculations have been made for the 
plasma density <n> = 0.7 x 10°° m® and the ablation 
region width A = 0.9. When Z = 1.2 the fusion power 
with charge 


stationary maintenance decreases 


approximately by 1 MW compared to the case with 
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Fig. 10 Radial distributions of the specific powers of 
bremsstrahlung gq, and cyclotron q, radiations in the case 
of fuel injection into the region of the plasma column 
radius middle part (x; = 0.5, Ph; = Phe = 5 MW, A = 0.24, <n> 
= 7x10" m®%, t= 36.2 s). 


MW 


T . T td T ° T " 1 


0. 5 10 15 2 25 30 35 40 45 
Time, s 
Fig. 11 DT-fusion power in the plasma heating process as a 
function of the time with different values of the Z plasma 
(A = 0.9, <n> = 0.7 x 10” m®, P); = Phe = 5 MW, Z = 1, 1.2, 
1.5). 
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Fig. 12 Plasma component temperature as a function of 
the time in the plasma heating process with different values 
of the Z plasma (A = 0.9, <n> = 0.7 x 107” m°, Py = Phe = 
5 MW). 

Z = 1. When Z = 1.5 the charge development leads to 
the fusion power increase, approximately to 20th 
second with a maximum value near 8 MW, followed by 
the smooth discharge decay. As Z = 1.5 the plasma 
temperature also is decreasing with time (Fig. 12). 
Such a plasma behavior is, probably, caused by the 
density profile realignment when the total loss power, 
with taking into account the bremsstrahlung power 
increase, exceeds the heating source power. At the 
plasma density 0.8 x 10°? m® the discharge decay 
begins just at Z = 1.2. 
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4. Discussion and Conclusions 


The possibility of neutron source development has 
been studied applying the stellarator system with linear 
dimensions increased by a factor of 1.5 compared to the 
experimental facility LHD. The calculations carried 
out under the assumption of the neoclassical transport 
show that using the 10 MW plasma heating sources, it 
is possible to obtain plasma parameters with which in 
the system the nuclear fusion power from one to 
several tens MW is produced. 

According to the investigation results the following 
can be concluded: 

(1) Optimum results, from the point of view of 
obtaining a maximum fusion power output, are reached 
by the fuel injection into the plasma column center. 
Besides, a more narrow fuel pellet evaporation region, 
a higher fusion power. 

(2) The fuel injection into the plasma center leads to 
the formation of strongly peaked plasma density 
distributions. 

(3) The nuclear fusion power increases with average 
plasma density value increasing. 

(4) The fuel injection into the peripheral region of 
the plasma column results in the plasma density spatial 
distributions, the maxima of which take place in the 
fuel injection region. 

(5) The impurity presence in the plasma causes the 
discharge decay just at Z= 1.5. 

In addition, it should be noted that the formation of a 
strongly peaked plasma density distribution is 
accompanied by the gas-kinetic plasma pressure 
increase when the value of / reaches more than ten 
percent that can create the problem for the stability of 
such distributions. 

The study has been carried out under the assumption 
of the neoclassical transport. The influence of some 
plasma loss mechanisms was not taken into account 
here. Therefore the question of calculation data 
application to the real system depends on the existence 
of another, not considered here, transport mechanisms. 

Furthermore, an additional study is necessary for the 


problem of breeding blanket location in the system 


under consideration. 
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Abstract: Orthogonal frequency division multiplexing (OFDM) is a special form of multi-carrier transmission that uses the policy of 
divide and rule. In this scheme, a large number of orthogonal, overlapping, narrow band sub-channels (subcarriers) are transmitted in 


parallel and divide the available transmission bandwidth. This techniqueis originally based on the Fast Fourier Transform of the 


information data. In order to improve the performance of the OFDM and overcome some limitations, an alternative OFDM approach 


based on the Wavelet Transform is proposed. In this paper, we study the performance of such systems in additive white Gaussian 
channel (AWGN). MATLAB simulations are realized and performance comparisons are presented. 


Key words: Multi-carrier modulation, discrete wavelet transform, wavelet packet transform, Fourier-based OFDM, wavelet-based 


OFDM. 


1. Introduction 


In wireless communications field, there is a rising 
demand for high performance, high capacity and high 
bit rate. Digitally techniques were adopted for these 
purposes. OFDM is a digital modulation technique that 
allows transmitting high data rates over extremely 
hostile channels at a comparable low complexity. In 
recent years, it has been adopted in a wireless local and 
metropolitan area networks and has been used in most 
of popular standards such as wireless fidelity (WiFi), 
Worldwide Interoperability for microwave access 
(WiMax), Digital Video Broadcast (DVB), 3rd 
Generation Partnership Project/Long Time Evolution 
(3GPP/LTE) standards. Due to worse transmission 
channel conditions, the interferences could attack the 
transmitted signal. The OFDM technique emerged to 
overcome this problem and combat channel distortion. 
It can be viewed as using many slowly modulated 
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narrowband signals rather than one rapidly modulated 
wideband signal [1]. It is a promising multi-carrier 
transmission technique, which divides the available 
spectrum into many subcarriers, each one being 
modulated by a low rate data stream, and then 
transmitted in parallel [2]. This division of the total 
bandwidth is an applicable way to deal with 
interferences that could attack the transmitted signal 
and to improve the spectral efficiency. Conventional 
OFDM is implemented using Fast Fourier Transform 
(FFT) and its counterpart Inverse FFT (IFFT) offering 
subcarriers orthogonality which prevents any expected 
carriers interference. FFT has a major drawback 
representing by high level side lobes. This leads to high 
interference and lower performance levels [3]. Thus, a 
guard interval or specifically cyclic prefix (CP) was 
inserted to mitigate this type of interference and to 
overcome the delay spread of the channel [1]. However, 
the CP brings a disadvantage of reducing the spectral 
efficiency for the transmitted signal by consuming 
bandwidth and power [4] without carrying additional 
information. Therefore, FFT based OFDM system 


The Performance Analysis of the Wavelet-OFDM New Scheme in AWGN Channel 101 


presents another major problem presented by the high 
peak to average power ratio (PAPR). 

An alternative OFDM design was proposed to 
overcome the limitations of the FFT based OFDM 
system. This is done by implementing a system with a 
discrete wavelet transform (DWT) instead of the FFT. 


The use of wavelet promises to reduce the 
inter-symbol-interference (ISI) and the 
inter-carrier-interference (ICI). It offers a high 


suppression of sides lobes. The Discrete Wavelet 
Transform is used in a variety of signal processing 
applications, such as video compression, object 
recognition and numerical analysis. Recently, the 
wavelet OFDM approach has been introduced into the 
communication field and it is included in the IEEE 
1901-2010 standard. Introducing of the DWT with the 
OFDM is extended to another form of wavelet 
transform; it is the Wavelet Packet Transform (WPT). 
Many researches have been carried out on this domain. 
The performance of Discrete Fourier Transform 
DFT-OFDM and DWT-OFDM using Binary Phase 
Shift Keying (BPSK) modulation are compared in 
Ref. [5] on three different channel models. While the 
impact of constellation size of Quadrature Amplitude 
Modulation (QAM) over FFI-OFDM and 
WPT-OFDM in AWGN channel is given in Ref. [6]. 
The performance comparison of DFT-OFDM and 
DWT-OFDM using arbitrary FIR channel with different 
types of wavelet packets is presented in Ref. [7]. In this 
work, we investigate the performance of wavelet 
OFDM system using the Haar wavelet, due to its some 
motivating features, upon an AWGN channel. 

The paper sections are organized as follows. Section 
2 presents the conventional Fourier-based OFDM. 
Section 3 theoretically describes the Wavelet-based 
OFDM. Simulation results are discussed in section 4. 
Finally, we conclude and present future work in 


Section 5. 
2. Fourier-Based OFDM 


In the FFT-based OFDM system, the information bit 


stream is first mapped using common digital 
modulation schemes such as Phase Shift Keying (PSK) 
and QAM schemes at the transmitter side. In this 
process, the serial binary data are converted into 
complex symbols representing constellation points. 
The Fast Fourier Transform (FFT) is applied at the 
receiver side to decode the signal. Then, each N 
symbols will be transmitted over N_ orthogonal 
subcarriers through the IFFT. The output of the IFFT 
block represents the OFDM symbol in the continuous 
time domain that is stated below: 


1 N-1 er 
(hey Xe" “04feT-- C1) 
In 2 - 


where, Af is the sub-carriers separation, XK is the FFT 
symbol information. 

After the OFDM symbol being constructed, a guard 
interval is added to mitigate the effect of ISI, which is 
caused upon propagating through the channel. Thus, a 
Cyclic Prefix is inserted at the end of each OFDM 
symbol (Fig. 1). At the receiver side, the data are 
reversely reconstructed using the FFT block in order to 
get the initial transmitted bit stream after applying the 
corresponding demodulation scheme. However, adding 
a CP reduces the spectral efficiency. Moreover, the 
high level side lobes create by the rectangular window 
of the FFT may cause interference in the frequency 
domain. These lacks have to be overcome. 


3. Wavelet-Based OFDM 


Unlike Fourier analysis which divides the whole 


bandwidth into orthogonal and overlapping subbands 
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of equal bandwidths, Wavelet transform (WT) is a 
multiplexing transmission method in which data is 
assigned to subbands having different time and 
frequencyresolutions. This allows a form of 
time-frequency analysis (or translation-scale in wavelet 
speaks). WT uses local orthogonal waves as_ basis 
functions to process data instead of sinusoids orthogonal 


basis. 
3.1 Discrete Wavelet Transform (DWT) 


The set of scaled and shifted wavelets constitutes an 
orthogonal basis of DWT, presenting an applicable 


characteristic invested in the subcarriers overlapping: 


Y mn(t) = 2? y(2"t—n) (2) 


where, y(t) denotes the mother wavelet, m and n are 
the scale and translate parameters, respectively. We 
define another function g(?) called the scaling function 
or the father wavelet where the shifting and scaling of 
this function will cover the entire square integral real 
space. y(t) and @(f) are used Q(t) together to process 
any signal and analyze it. They are related in the 
following relations: 


(t) = > h,(n)V29(2t—n) (3) 
w(t) = >°h, (n)W2g(2t-n) (4) 


hg(n) and hy (n) are respectively a low pass and a high 
pass filter banks. The low pass filter gives approximate 
information about the signal while the high pass filter 
produces detailed information. We mainly interest in 
the approximate coefficients where we can find most of 
the information. Fig. 2 shows the process done with 
DWT filter banks. The low or approximate coefficients 
are down sampled and re-filtered. 

The signal reconstruction is done in the same way by 
filtering the coefficients with the same filter banks used 
for decomposition-Note that both filters are also 
assumed to have perfect reconstruction property, this 
means that the input and output of the two filters are 
expected to be the same [8]. Finally, they are 


up-sampled to get the original signal. 
3.2 Wavelet Packet Transform (WPT) 


In the DWT, only the approximate coefficients are 
decomposed again and decimated by 2. WPT both 
approximate and detailed coefficients undergo the 
same operation of decomposition and decimation at 
each level [3]. For the wavelet packet, we can decide 
how to decompose the high and low frequency parts, 
after each decomposition, we can decide whether to 
decompose the signal in the low/high frequency 
domain or not. They form basis, which retain many of 
the orthogonality, 
localization properties of their parent wavelets. 


smoothness, flexibility, and 


3.3 Wavelet Transform-Based OFDM 


Wavelet localization feature can provide much lower 
side lobes magnitude than those of Fourier transform 
(FT) and the loss of orthogonality leads to lesser ISI 
and ICI. 
frequency synchronization to get a good performance. 


Therefore, OFDM requires time and 
This can be overcome by the use of the WT-based 
OFDM approach. 

The FFT block in the original diagram is now 
replaced by the DWT block. The symbols resulted 
from constellation mapping undergo inverse discrete 
wavelet transform (IDWT) operations and construct 
the desired OFDM signal. The operation is done by the 
aid of quadrature mirror filters that are the low pass 


Low Pass 


¥ 
High Pass 


High Pass 


Fig. 2 Schematic diagram of a discrete wavelet transform. 
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IDWT filter banks. Perfect 
reconstruction and orthonormal basis properties are 


and high pass 


satisfied, where symbols are convoluted with those 


filters and then multiplexed together forming the signal. 


The data are pointed as the approximate coefficients 
of the filters and the detailed coefficients are 
zeropadded [8]. At the receiver where the DWT block 
replaces the FFT one, the signal is_ reversely 
decomposed producing the approximate and_ the 
detailed coefficients. Actually, the useful information 
is contained in the approximate coefficients. QAM 
demodulation is then done and the data stream is 
retained. 

Similar to IFFT, IDWT matrix is a transpose of the 
original. Both transforms can be viewed as a rotation in 
function space to a different domain [5]. For the FFT, 
this new domain contains basis functions that are sines 
and cosines. For the wavelet transform, this new 
domain contains more complicated basis functions 
called wavelets, mother wavelets, or analyzing 
wavelets. 

Thus, an improved performance must occurwhen the 
wavelet transform is introduced within — the 
multi-carrier OFDM-based system. However, the 
logarithmic division of the bandwidth in the wavelet 
transform is not well suited for OFDM multi-carrier 
communication [9]. This motivates us to introduce the 
WPT. In addition, WPT is more robust to the tone and 
impulse interference added into a transmission packet. 
So the transmission throughput is kept in WPT even in 
these severe interference environments [10]. However, 
while the DWT is less computationally complex than 
the FFT due to the logarithmic division of frequency, 
WPT is more complicated than FFT. This is the price 
that must be paid. 


4. Simulations Results and Discussions 


Theoretically, WT promises an improved 
performance in comparison with the FFT. This fact has 
to be exploited for the OFDM systems. However, this 


idea is still incompletely studied and it needs to be 


proved in practice. Thus, we are motivated to compare 
the performance of the FFT-based OFDM system with 
the DWT-based OFDM one. The MATLAB 
simulations were done to achieve this goal. 

Firstly, the conventional FFT-based OFDM block 
diagram is designed using 128 subcarriers. Then the 
Discrete Fourier Transform (DFT) is replaced by the 
DWT. As the Mother Wavelet, the HAAR wavelet is 
used. The input data stream is modulated with the 
16-QAM (Quadrature Amplitude Modulation). The 
channel is assumed as Additive White Gaussian Noise 
(AWGN) channel. The Bit Error Rate (BER) and the 
Power Spectral Density (PSD) are the main keys used 
to comparison. It is important to note that the choice of 
the wavelet transform is very critical since it is related 
to the orthogonality property. For some wavelet 
functions, bases are not orthogonal to each other. This 
causes the AWGN to be correlated within a subcarrier 
and an AWGN channel does not remain an AWGN one. 
Thus, the result is worse than the single carrier case. 
For the Haar wavelet, it is not the case and no timing 
frequency synchronization is done in Ref. [9]. 
Moreover, it is more adapted to reduce the PAPR [10]. 


4.1 Bit Error Rate Performance Criterion 


The BERs of the OFDM system transmission are 
calculated over a several Monte Carlo runs and are 
plotted with respect to the signal to noise ratio (SNR) 
or E,/No (dB). Table 1 represents the parameter values 
used for simulation. 

Fig. 3 shows the performance comparison between 
FFT and DWT OFDM system. We can clearly 
observethat the DWT based system outperforms the 
FFT one and provides a “3 dB” performance gain. 


Table 1 System parameters. 


; FFT-OFDM DWT-OFDM 
Variable 2 : 
Matrix value Matrix value 

Numberofsubcarriers 128 128 
CP 26 0 
Bitstream 512 512 
Repetitiontimes 20,000 20,000 
Waveletused - Haar 
SNRrange(dB) 1-20 1-20 
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Fig. 3 FFT-based OFDM vs DWT-based OFDM. 


In order to observe the effect of the OFDM system, 
we compare the raw 16-QAM system to the 16-QAM 
OFDM system. First, using the FFT tool (Fig. 4), and 
then using the DWT (Fig. 5). The results show that the 
16-QAM DWT-based OFDM system outperforms the 
raw 16-QAM one with a 3dB as a gain, unlike in the 
case of the FFT-based OFDM system where the 
performance is too close to the raw-QAM system. We 
must note here that the main cause of this observation 
is that we assume the AWGN channel case which does 
not reflect clearly the OFDM benefits. But, we focus in 
this work on the effect of the WT with comparison to 
FT within the OFDM systems. 

Therefore, we can say that either QAM alone or 
QAM with FFT-based OFDM will do the same work 
under same conditions and both are less applicable 
when compared to DWT based OFDM system. 


4.2 Power Spectral Density Comparison 


Now, we compare in the frequency domain. The 
PSD for each transmitted signal is presented in Fig. 6. 
QAM is theoretically a single carrier modulation 
scheme, where the signal is transmitted over one single 
frequency. Thus, the spectrum of the 16-QAM should be 
a single carrier signal as shown by simulation (Fig. 6). 
Merging QAM with OFDM creates a multi-carrier 
transmission. Indeed, the FFT-based OFDM spectrum 
proves that, as shown in Fig. 6 where the variation of 
the power levels indicates the separation between the 
sub-channels. 


[| “EB 16-040 System 
| | =P 16QAM FFT-OFOM 
t—=—=—= 


12 snr (dp) ‘4 
Fig.5 16-QAM and 16-QAM with DWT-based OFDM. 


Now, when substituting the FFT with DWT, as we 
proposed in this project, we should maintain the 
OFDM ulti-carrier principle. This is theoretically 
achieved and the studies about Wavelet Transform 
confirm this concept. Mathematically, we can prove 
that the use of IDWT leads to a multi-carrier 
transmission. 

The operation undergoes a set of translations and 
scales of the wavelet function along the whole signal 
and convolving the scaled and shifted wavelet with our 
concerned signal. This leads to a set of subcarriers 
orthogonally spaced [11]. As a result, the PSD for both 
FFT-based system and DWT based one should 
logically have the same form which meets with the 
OFDM theoretical one. Well, this is the main concept, 
but does this really happen with the MatLab 
simulations? 

Unfortunately, simulations do not show a 
coincidence of the DWT-based OFDM spectrum 
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Fig.6 Power spectral density of (a) QAM; (b) FFT-based 
OFDM and (c) DWT-based OFDM system. 


(Fig. 6c) with the FFT-based one (Fig. 6b). Therefore, 
this PSD seems different from the theoretical spectrum 
and appears as a single carrier transmission. Otherwise, 
the figure shows a gain in the spectral efficiency 
offered by introducing of DWT: the PSD width of the 
DWT-based OFDM is less than that of the FFT-based 
OFDM since they are not adding a CP. 

In fact, with Mat Lab “DWT”, the convolution 
between the wavelet function and the data symbols is 
done only one time, nor shifting neither scaling are 
applied on the wavelet. This means that the principle of 
multi-carrier modulation is not applied since the reason 


behind the creation of subcarriers is not verified. 


WPT-OFDM psd 


a a 
VALS FTA 

a 

| mip 


N 
= 
oO 
= 
> 
rs) 
e 
o 
5 
& 
£ 
é 
a 


1.5 
Frequency (kHz) 


Power spectral density of WPT based OFDM 


Fig. 7 
system. 


PETE] 
| @= FFT-OFOM 
| “8 weT-OFoM 


Fig. 8 WPT-based OFDM and FFT-based OFDM BER 
performance. 


Thus, MatLab is not the right way to get the exact 
answer for the comparison between FFT-based OFDM 
and DWT-based OFDM since the first one is a 
multi-carrier modulation type and the second one is a 
single carrier. Raw QAM and DWT-OFDM can be 
compared; they are both a single carrier modulation 


types. 
4.3 Wavelet Packet Transform (WPT) Performance 


Using the 'DWT' function does not lead to a 
multi-carrier OFDM system. An alternative way can 
overcome this lack: this is the WPT. Analytically, the 
WPT is the Real Discrete Wavelet Transform (DWT). 
This MATLAB function leads to multi-carrier OFDM 
system. Thus, we must obtain the same PSD of the 
traditional OFDM system. The simulation proves this 
theoretical fact, as shown in Fig. 7. To reach a complete 
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answer about this proposed solution we calculate the 
BER under the same conditions and we obtain the Fig. 
8. As it is clear, the WPT has the same performance of 
the FFT-OFDM system and does not provide any gain. 


5. Conclusions and Future Work 


The 
OFDM.-transmission performance by introducing the 
DWT mathematical tool instead of the FFT one. 


Indeed, the Wavelet Transform outperforms the 


aim of this work is improving the 


Fourier transform in communication 


applications. MATLAB simulations show that this 


many 


idea provides “3 dB” in term of BER gain. However, 
the system lost his multi-carrier characteristic. This 
observation motivates us to introduce the Wavelet 
Packet Transform WPT as an alternative solution. 
This way leads to a multi-carrier OFDM system, but 
does not present any performance gain by comparing 
to the traditional OFDM system. This result leads us 
to a conclusion that the comparison of FFT-based 
OFDM and WT-based OFDM using MATLAB is not 
the right way to get the exact answer for our proposal. 
It is better wise to compare raw QAM with the DWT 
system; they are both a single carrier modulation. 

We propose to perform this study using another 
simulation way to meet results with the theoretical 
properties. As practical case, it is suitable to study the 
case of the fading transmission channel. 
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Abstract: Biofilm is a community of bacteria, less susceptible to traditions treatments. Although photodynamic therapy (PDT) is a 
very effective way to microorganism inactivation, in biofilm it is not as efficient as it is in planktonic bacteria cultures. The increment 
of an element to increase the effectiveness of PDT was our aim. Therefore, this in vitro study evaluates the susceptibility of a biofilm 
formed by Streptococcus mutans on metallic surface of orthodontic accessories under the application of PDT with a surfactant. Samples 
obtained from blades of orthodontic bands (NiCr), where used as adhesion surface for the biofilm. They were treated with 1 mg/ml of 
curcumin, with 0.1% of sodium dodecyl! sulfate and exposed to 30 J/cm? of light (455 nm). Eight experimental groups were studied, 
including the positive and negative controls. The results show that the group with PDT and surfactant had a significant decrease (p < 
0.001) in viability. In this case, the reduction observed was of 5.6log10 (CFU/ml) in comparison to the control group. We have shown 
that, even though the biofilm is very tough and complex structure, we are able to promote almost the complete inactivation of S. mutans 
in systems similar to an orthodontic treated patient’s mouth. 


Key words: Photodynamic therapy, light-induced damage, oral biofilm, streptococcus mutans, curcumin, surfactant. 


1. Introduction For that reason, biofilms are precursor for most 


; a ; common oral diseases, such as caries and periodontal 
There are two forms for bacteria to survive in their ; : : ; 
: : disease. Prior to the formation of the biofilm, a thin 
natural environment. They may be free-floating, as a : : 
; . on microorganism layer adheres on the surface of the tooth 
planktonic cell, or attached to a surface, as in a biofilm. ; : : ; 
ae , ; or gum. These microorganisms use salivary proteins as 
The biofilm is a complex community in which most of ; 
= ; . ; substrate for adhesion. Streptococcus mutans produces 
bacteria live. It is a well-organized community that ; ; : 
; ; adhesins that together with these proteins are used to 
adheres to surfaces and is embedded in an . 
; : form the substrate of adhesion. The extracellular 
extracellular adherent layer. Bacterial resistance ; 
: . — substrates that they produce allow the accession of 
increases when they are organized as biofilm. Also, ; ; 
; . cae j other types of microorganism to form the mature 
once in a biofilm, the bacteria display different as ; 
_ biofilm, and finally the dental plaque. 
characteristics from those that they had as a ae ; : 
. : — . The dental biofilm can be categorized into sub- and 
free-floating organism. Biofilms can be found in ae ; j 
supra-gingival. It is formed subsequent to an increase 
many places, and the human mouth has perfect : ; ee : 
S ; of saccharolytic and acidogenic microorganisms, such 
conditions for the bacteria to grow. ; ; 
as Streptococcus mutans and lactobacillus, which 
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demineralization and caries. Four components are 
needed to create the carious lesion [1, 2]: a host, 
high-carbohydrate diet, microorganisms and time. 

The Streptococcus mutans, that may be one of the 
most common bacteria in the oral cavity, produces a 
polyhedral matrix to structure the biofilm. The 
microorganisms existent in a biofilm become sheltered 
by this polyhedral matrix, reducing the action of 
antimicrobial agents. Compared to planktonic cells, the 
biofilm structure protects the bacteria, and other 
microorganisms. The thickness of this structure 
prevents contact between antimicrobial agents and 
microorganisms, keeping the chemicals away from the 
deeper layers of the biofilm. As a result, the biofilm 
boosts the growth rate of the microorganisms [2-4]. 
Consequently, microorganisms embedded in biofilms 
are 10 to 1000 times more resistant to antimicrobial 
agents [5, 6]. 

One cause that enhances the biofilm formation in the 
mouth is the installation of orthodontic accessories [7]. 
The usage of orthodontic accessories creates new 
retention surfaces for microorganisms in the oral cavity, 
increasing their retention. Therefore, it increases 
biofilm formation and the amount of bacteria in the 
mouth [7]. The high level of S. mutans in the oral cavity 
of patients with orthodontic accessories, compared to 
non-orthodontic treatment patients, is cause both by the 
presence of the accessories and also by the poor 
hygiene control of it by the patient. 

Photodynamic therapy (PDT) is a well-known 
medical therapy for cell inactivation and 
microorganism control [8, 9]. This technique is based 
on the combined use of light and some types of drugs, 
called photosensitizers (PS), which produce an 
oxidation reaction. This oxidative reaction can them 
cause the death or inactivation of selected cells. PDT 
can also be used as a therapy with a potent 
effect [10]. 


microorganisms by PDT can ensure several advantages 


antimicrobial Inactivation of 


over the traditional use of antimicrobials. First, 


bacterial killing is rapid, reducing the need to maintain 


high concentrations of chemicals for long periods, as in 
the case of antibiotics and antiseptics usage [11]. 
Secondly, bacterial killing is not linked to the 
intervention of chemicals, so, development of 
resistance is improbable to happen [12]. Finally, since 
neither the PS, nor the light, are bactericides by 
themselves, the destruction of bacteria can be 


controlled by restricting the irradiated region, 
preventing the destruction of the normal micro-biota 
elsewhere. 

PDT is an alternative therapy for preventing and 
treating dental caries and periodontal diseases. It is an 
innovative way to control the formation of the bacterial 
biofilm, controlling the incidence of these pathogens. 
Most studies present high efficiency of PDT to 
planktonic bacteria. On the other hand, since the 
biofilm protects its microorganisms, the results over 
clinical trials are not as efficient [13]. Therefore, 
researches over systems that emulate the patients’ 
mouth are so important to find the “perfect therapy”. 
Thus, the purpose of this paper is to evaluate the effect 
of the antimicrobial activity of PDT on biofilms formed 
on orthodontic accessories. This in vitro study 
evaluates the use of PDT and a surfactant on metal 
surfaces, observing the susceptibility of biofilms 
formed by S. mutans. Our objective is to increase the 
performance of PDT, adding one element to the 
treatment, in a system that emulates the mouth real 


environment. 
2. Material and Method 


For this study 24 metal samples measuring 5 x5 
mm were used, obtained from orthodontic band blades 
(NiCr). On the metal surface, S. mutans biofilms were 
growth and threated with PDT and a surfactant. The 
PDT was induced by curcumin and light emitting 
diode (LED) in the blue range, and the surfactant used 
was the sodium dodecyl sulfate (SDS). There were 
tested eight experimental groups, each test 
repeated three times and each solution measured 


twice. 
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2.1 Biofilm Growth and CFU Preparation 


For Streptococcus mutans biofilm, an inoculum 
from stock culture (ATCC 25175) was grown in Brain 
Heart Infusion (BHI) broth (Difco, Detroit, USA) and 
incubated at 37 °C for 24 h. The suspensions were 
centrifuged, 3,000 rpm per 15 min (Excelsa II 
centrifuge, model 206-BL, FANEM), and the bacterial 
pellet was dispersed into BHI broth with 20% sucrose. 
To archive the desired population density, the culture 
was adjusted to obtain standardized suspension 
containing 10° cells/ml, the optical density of the final 
suspensions was _ verified 
spectrophotometer (FEMTO 600). 

The biofilm was grown in orthodontic band metal 


using a digital 


blades (orthodontic appliance, Tecnident, Sao Carlos, 
Brazil), cut into 5 x 5 mm squares, and sterilized in 
autoclave. Each sample was placed inside individuals 
wells of a 24-well microtiter plate. Then, aliquots of 
1 mL of the S. mutans were transferred to the plates. 
The biofilms were incubated at 37 °C up to 7 days, with 
the growth medium changed every two days. All the 
samples were washed in PBS solution and placed in a 
new well into the 24 microtiter plate to be prepared for 
experimental treatment. They were separated between 
eighth groups and treated accordingly. 

After treatment, the samples and theirs solutions 
were transferred to eppendorfs, 1ml of PBS was added, 
and they were homogenized, to break and loosen the 
biofilm into the solution, during 8 min in an 
ultrasound shaker (150 Watts, 50 Hz), before the 
metal squares were removed. The solutions were 
diluted up to 10° and uniformly spread to petri dishes 
containing culture medium BHI AGAR (Difco, 
Detroit, USA) plus 20% 
aerobically cultured at 37 °C for 48 h prior to the 


sucrose. They were 


colony-forming unit (CFU) count. 
2.2 Experimental Treatment Procedure 


The experimental groups included: control group 
(P-L-S-); light (P-L+S-), surfactant (P-L-S+) and 
(P+L-S-) groups; 


photosensitizer surfactant-light 


(P-L+S+), photosensitizer-surfactant (P+L-S+) and 
PDT (P+L+S-) groups; and the PDT-surfactant 
(P+L+S+) group. Each group was named accordingly 
to the addiction of the photosensitizer (P+ or P-), 
surfactant (S+ or S-) and application of light (L+ or L-). 

All the samples were placed in a new well into the 24 
microtiter plate and received the correct solution 
necessary for treatment. To the control (P-L-S-) and 
light (P-L+S-) groups, there was added 1 ml of a 
solution of PBS with 5% of DMSO. There was added 
1 ml of surfactant solution to the surfactant (P-L-S+) 
and surfactant-light (P-L+S+) groups. The surfactant 
solution was made in PSB with 0.1% of SDS and 5% of 
DMSO. To the photosensitizer (P+L-S-) and PDT 
added 1 mL of a 
photosensitizer solution. This solution had 1 mg/ml of 
curcumin, and was made in PBS with 5% of DMSO. 
And to the photosensitizer-surfactant (P+L-S+) and 
PDT-surfactant (P+L+S+) groups, there was added 1 
mL of a photosensitizer-surfactant solution, which had 


(P+L+S-) groups it was 


1 mg/mL of curcumin, 0.1% of surfactant, and 5% of 
DMSO in PBS. 

All the experimental groups were kept in the dark for 
5 min, incubation of the PS, before the application of 
light in the designated groups. The experimental 
groups with light (L+) were then irradiated in a 
home-made blue LED device (high power royal blue 
LEDs, 355mW, centered at 455 nm with 20 nm 
bandwidth), for 10 min and 54 s, ensuring a 30 J/cm? 
dose at all the wells. The LED device (Fig. 1), called 
Biotable (LAT-IFSC-USP), was made to guarantee 
uniform distribution of light (46 mW/cm? intensity) in 
the plane were the samples were placed. The 
experimental groups (L-) where no light was applied 
were kept in the dark at room temperature, for the same 
amount of time. 


2.3 Data Analysis 


The effect of the photodynamic therapy and 
surfactant on the biofilms was evaluated by counting 
colony-forming units. The results were statistically 
evaluated by ANOVA, differences were considered 
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Fig. 1 Photography of the 24-well blue Biotable and details of the metal plate samples, placed in a 24-well microtiter plate, 


being put in place on the LED device. 


when p < 0.05. The difference between treated groups 
with the control group (P-L-S-) was made with the 
OneWay ANOVA with post hoc Turkey. 


3. Results and Discussion 


To ensure biofilm formation, selected metallic 


samples were imaged with confocal microscopy (Fig. 2). 


The images showed the formation of S. mutans biofilm 
after 7 days of cultivation on the metallic samples of 
orthodontic appliance (NiCr). The photodynamic 
inactivation was also observed, by means of a live-dead 
biological marker, imaging plate before and after PDT 
treatment. The confocal reflection images show the 
metal surface, indicating the presence of colonies of S. 
mutans. The confocal fluorescence images show the 
fluorescence of the propidium iodide (PI) dye and 


curcumin in the bacteria cells. The colonies can be 
highlighted when both images are superposed, showing 
the correspondence between the colonies seen by the 
reflection and fluorescence images. 

The confocal reflectance images present the surface 
of the metal sample. It was expected that the 
orthodontic appliance squares show some roughness, 
and appears in the image as a random pattern. It is also 
noticed a pattern similar to a bunch of grapes, and the 
superposition of the reflectance and fluorescence 
images confirm that those patterns correspond to the 
biofilm. Meanwhile, the fluorescence image show red 
and yellow-green dots. The yellow-green dots are the 
fluorescence of the curcumin within the bacteria cells. 
The red dots are the fluorescence if the propidium 


iodide (PI) dye, which appear in the dead cells. 


Fig. 2 Biofilms of S. mutans on metallic sample with 7 days, after PDT treatment: (a) Confocal reflection imaging, the image 


show the appearance of the metal surface, indicating the presence of S. mutans colonies; (b) confocal fluorescence image, 
showing the fluorescence of the propidium iodide (PI) dye and curcumin. The colonies appear as bright spots in the image, red 
indicates dead cells and yellowish-green show the living microorganisms and (c) reflection and fluorescence images superposed 
(a + b). The image shows the correspondence among the colonies seen by reflection and fluorescence images. The scale bar 


corresponds to 20 pm. 
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The bacterial effects of PDT and the surfactant were 
evaluated by the cell viability (CFU/ml). The efficacy 
of the treatment differs depending on the association of 
curcumin and surfactant (Fig. 3). They showed 
significant reductions in the viability of the S. mutans 
(over 1 logi0) for all groups, except for the 
experimental conditions (P-L-S+) and (P+L-S-). The 
highest reduction in the cell viability of the biofilms 
was equivalent to 5.6 log10 (CFU/ml) and correspond 
to the association of PDT and surfactant. 

The results show remarkable behavior of the bacteria 
population depending on the association of curcumin, 
surfactant and the application of light (Table 1). In some 
cases the illumination can promote biostimulation, but 
none of the studied cases presented a relatively relevant 
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increase on the bacterial population. Applying just PDT 
(P+S-L+) or surfactant (P-S+L-), it was observed small 
decreases on the S. mutans population. The same occur 
when light was applied mutually with surfactant 
(P-S+L+), but in this case the antibacterial effect was 
somehow smaller. The test group where only the 
photosensitizer was used (P+S-L-) was statistically 
similar to the control group, confirming that curcumin, 
by itself, is not capable to a significant bacteria 
inactivation. On the other hand, the association of 
curcumin and SDS (P+S+L-) upheld an even smaller 
decrease in the population. However, when the three 
where used together (P+S+L+), it was observed a huge 
antibacterial effect (-5.6 log10). In this case, only one in 
a thousand part of the population survived. 
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Fig.3 Reduction of CFU for S. mutans biofilms. Individual difference made by ANOVA One Way: post hoc Turkey. Groups 


with the same letter are not significantly different (p > 0.005). 


Table 1 Expression of cell viability (CFU/ml), related to the control group data, for all experimental groups. Groups with the 


same letter are not significantly different (p > 0.005). 


Experimental groups CFU percentage difference Expression (log10) 
(P-S-L-) * 0 

(P-S-L+) * +29.17 +0.1 

(P-S+L-)? -99.22 -2.1 

(P-S+L+) °° -97.84 -1.6 

(P+S-L-)*4 -42.91 -0.2 

(P+S-L+) °° -97.46 -1.5 

(P+S+L-) °7° -93.37 -1.1 

(P+S+L+) ! -99.999792 -5.6 
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The oral microbiota is composed of more than 500 
types of microorganisms, even though S. mutans is the 
etiologic agent of caries and demineralization. It is also 
responsible for the adhesion, colonization and dental 
biofilm formation [14]. This 
particularly the susceptibility of S. mutans biofilm on 
the metal surface (NiCr) of orthodontic blade band by 
of photodynamic therapy with curcumin 
(1 mg/mL), 0.1% of SDS and LED illumination. As 
expected, we have shown that photodynamic therapy 


study evaluated 


means 


reduces the number of living cells of S. mutans in a 
biofilm model created on metallic surface [15-17]. The 
presence of photosensitizer, without irradiation, was 
not able to achieve an effective antibacterial outcome. 
This confirms the essential rule of light to obtain an 
anti-microbial effect [18, 19] by means of PDT. 

The most significant finding was the lack of viability 
of S. mutans on the surfaces of the samples treated with 
curcumin combining PDT with surfactant. The result 
for this group (P+S+L+) is above the minimal 
bactericidal concentration, which is the minimum 
concentration able to reduce 4 log10 on the survival of 
a microorganism [20]. Significant decrease on the 
viability of microorganisms was also observed when 
biofilms were exposed to the surfactant alone (P-S+L-), 
surfactant with light (P-S+L+), photosensitizer with 
light (P+S-L+), or the combined use of surfactant and 
curcumim without light (P+S+L-), compared to the 
control group (P-S-L-). For these groups, the reduction 
on the bacteria viability was above 90%. Depending on 
the application, 90% of reduction on a microorganism 
population might be sufficient. Although the results 
show that the bacterial inactivation can be almost 
complete when the surfactant was used together with 
PDT. 

The surfactants are amphipathic molecules 
composed of a hydrophobic portion and a hydrophilic 
portion. The first is often a nonpolar hydrocarbon chain, 
while the second can be ionic (cationic or anionic), 
nonionic or amphoteric. The surfactants, for their 
biochemical characteristics, are highly potent and act 


as emulsifiers, reducing interfacial and superficial 
tension [21]. They are responsible for modifications on 
the ion channels of the bacterial membrane via 
trans-membranes proteins, similar to the porins. This 
allows the passage of solutes through the membrane, 
increasing the permeability of these membranes [21-24]. 
Surfactants also work as: emulsifiers, dispersing in 
water hydrocarbons or other insoluble compounds; 
reducing the adhesion and releasing surfaces cell; and 
presenting antibiotic activity. 

The mechanism of action of PDT is not connected to 
the mediation of chemical radicals acting in a single 
target, which is the case of the antibacterial products. 
Photosensitizes act producing reactive oxygen species 
(ROS) which have no specific site of action. This 
avoids the formation of strains resistant to 
antibiotics [20, 21, 25, 26]. Photodynamic therapy 
applied as antimicrobial therapy can be _ also 
advantageous over mechanical removal methods, since 
it can reach places such as recesses and protrusions of 
the orthodontic accessories [27]. It may also enable 
safe treatment of patients with special needs or 
difficulty in oral hygiene during orthodontic treatment, 
avoiding oral infections, such as dental caries, 
gingivitis and periodontal pockets [13, 27]. 

Several hypotheses could explain the synergism of 
the PDT action with surfactant as antimicrobial therapy 
in S. mutans biofilm. The surfactant could be 
increasing the number of working cycles of PS, 
protecting it from oxidation, and increasing the 
efficacy of PDT. It could also be acting increasing the 
permeability of the curcumin through the cell 
membrane, increasing the concentration of PS 
molecules in the bacteria. One other hypothesis is that 
the surfactant is breaking the biofilm structure, and this 
allows the PDT to be more effective. Although the 
surfactant, SDS, can be used as antimicrobial by itself, 
the combined effects of it with PDT cannot explain 
such efficacy of the treatment. The reduction on 
bacteria viability was larger than one could expect from 


using both treatments separately. Therefore more 
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studies must be done in order to thoroughly understand 
this mechanism. 


4. Conclusions 


The result of this study shows higher reduction in the 


S. mutans population, compared to previous 
investigations [15, 19]. We demonstrate significant 
removal of S. mutans biofilm after application of PDT 
with SDS as antimicrobial therapy. PDT associated to 
surfactants promoted a high inactivation of S. mutans 
biofilm (5.6log10) on the surface of metallic orthodontic 
accessories. We have shown that, even though the 
biofilm is very tough and complex structure, we are able 
to promote almost the complete inactivation of S. 
mutans in a system similar to an orthodontic treated 
patient’s mouth. And this inactivation was achieved 
with a protocol that does not promote bacterial 
resistance that applied PDT and SDS. 

Given the findings, our results suggest that the 
surfactant acts on the microorganisms’ membrane 
making them permeable to solutes, and breaking the 
biofilm structure. The surfactant leads to higher the 
input of oxygen molecules and curcumin from the 
extracellular medium into the microorganism. 
Therefore, the synergistic action of the SDS and PS 
produces higher amounts of reactive oxygen species 
(ROS) upon irradiation, optimizing the antimicrobial 
effect of PDT. At the same time, protecting the PS 
molecules from oxidation and increasing its lifetime of 
the PDT cycle. These hypotheses should be put to test 


in supplementary investigations. 
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Abstract: In this work, enhancement of the light extraction efficiency of a 590 nm AlGalInP light-emitting diodes (LED) with a 
reflective top electrode (RTE) was investigated. A distributed Bragg reflector (DBR), consisting of AlAs/AlGaAs pairs, grown on an 


AlGaInP structure was used as a reflector for a reflective top electrode. It was found that a higher output power was observed from the 


AlGaInP LED with a RTE than from a conventional one. In addition, it was noted that the improvement in the output power depends 


strongly on the reflectivity of the reflector and that it exhibits a more effective performance with low injection currents. The increase in 


the optical output power was attributed to the enhanced extraction efficiency caused by a reduction of light absorbed from the emission 


region to top electrode through the RTE. 


Key words: AlGalnP, reflective top electrode (RTE), light-emitting diode (LED). 


1. Introduction 


III-V compound semiconductors are known as an 
attractive material for optical devices such as solar 
cells [1], laser diodes [2], and light-emitting diodes [3] 
amongst others. In particular, light-emitting diodes 
based on AlGaInP materials, with a wavelength 
ranging from 570 nm to 630 nm, have been developed 
for commercial applications such as illumination, 
traffic lights, vehicle lighting and indicators [4]. Ever 
higher efficiency AlGaInP LEDs have been required 
for high power devices with a low operating voltage. 
The efficiency of AlGaInP LEDs has been significantly 
improved by either increasing the internal quantum 
efficiency by inserting multiple hetero-structures and 
multiple quantum barriers, or by increasing the light 
extraction efficiency by using a distributed Bragg 
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reflector [5-7], surface texturing [8], bottom 
window [9] or nanoscale materials deposited on their 
surfaces [10]. 

Although numerous studies have been reported for 
AlGaInP LEDs with a high efficiency, certain intrinsic 
limitations still exist that considerably decrease their 
overall efficiency. In particular, it was shown that a 
significant reduction in the efficiency depends strongly 
[11, 12], as the 


light-emitting region in the active region is generated 


on an opaque top electrode 


mostly beneath the opaque top electrode. 

In this study, the fabrication and effects of a 
reflective top electrode on the light extraction 
efficiency of a 590 nm AlGaInP LED were 
investigated. This reflective top electrode was 
fabricated by inserting a AlAs/AlGaAs distributed 
Bragg reflector between the top electrode and the LED 
structure using a MOVD system. It was found that the 
light extraction efficiency of 590 nm AlGaInP LEDs 
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could be increased by using the reflective top electrode 
developed. 


2. Experiments 


The AlGaInP LEDs used were grown on an n-type 
GaAs (111) substrate at a 10° forward tilt [0-1-1] + 
0.5° using an MOCVD system. Here, trimethylgallium 
(TMGa), (TMAI), and 
trimethylindium (TMIn) were used as the group III 


trimethylaluminum 


sources; arsine (AsH3) and phosphine (PH3) were the 
group V sources; and disillane (SipH6) gas and 
cyclo-peantadienylmagnesium (Cp2Mg) were the n- 
and p-doping sources. Hydrogen (Hz) was used as the 
carrier gas for all sources. 

The cross-sectional schematics of the AlGaInP LED 
structures developed are shown in Fig. 1. In brief, the 
substrate is a GaAs lattice matched with both AlGaInP 
and AlGaAs. The AlAs/AlGaAs distribution Bragg 
reflector (DBR) is used to reflect upward the light 
emitted from the active region to the substrate. The n- 
and p-AlGaInP materials were used as the n- and 
p-confinement layers, respectively, by varying the 
component ratio of Al and Ga in the AlGaInP materials. 
The active region consists of 30 pairs of MQWs, 
having 10 nm thick AlGaInP wells and 30 nm thick 
AlGalInP barriers. In conventional AlGaInP LEDs, the 
GaAs top layer of 80 nm thick was grown on a p- 
confinement layer for ohmic contact. In the developed 
AlGaInP LEDs with a reflective top electrode, the 
AlAs/AlGaAs top distributed Bragg reflector 
(TDBR) , consisting of ~68nm thick AlAs and ~57nm 


@) (b) 
Fig. 1 Schematics of structures of (a) conventional 


AlGaInP LED and (b) AlGaInP LEDs with a reflective top 
electrode. 


thick AlGaAs periods, 


confinement layer for both reflection and ohmic 


was grown on the p- 


contact. For further improved information, the TDBRs 
of 4 and 12 periods were trialed for the reflective top 
electrode. 

In this work, all the epitaxial wafers used were 
divided into 9 mil chips (230 um x 230 um) using 
dicing system. 


3. Results and Discussion 


Fig. 2 presents the photoluminescence (PL) 
measurements for the AlGaInP LED epitaxial wafers 
with difference structures, which were obtained using a 
525 nm laser diode source, showing that all of these 
epitaxial wafers have the same wavelength of 590 nm. 
The figure shows that the intensity of the AlGaInP 
LEDs was decreased remarkably by the insertion of the 
TDBR, with the highest intensity being obtained in the 
AlGaInP LED epitaxial wafer without a TDBR. 
Additionally, a significant reduction in the intensity 
close to zero was observed from the AlGaInP LED 
with the TDBR of 12 periods. From these results, it 
was confirmed that the light emitted upward from the 
LED was reflected effectively through the use of a 
TDBR. 

A schematic of the structure of an AlGaInP LED 
with a reflective top electrode is shown in Fig. 3a. 
it can be assumed that 


From the schematic, 


the emission region in the active layerwas generated 


0.05 


—*— Conventional LED 
LED with TOBR of 4 penoas 


SLED with TOBR of 12 partods 
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Fig. 2 Photo-luminescence measurement of AlGaInP LED 


epi-wafer with different structures. 
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Fig. 3 (a) Light emitting path on lights emitted upwardly 


from emission area in AlGaInP LED chip with the RTE and 
(b) graphic of reflectivity and angle of incidence for the 
different DBRs used as reflector. 


mainly below the top electrode, and that most of the 
light emitted upward from the emission region was 
moved to the top electrode. Therefore, the efficiency of 
a conventional LED must be low due to the top 
electrode being absolutely opaque to the light injected. 
On the other hand, an even higher efficiency could be 
obtained from an AlGaInP LED with a reflective top 
electrode because light emitted from the emission area 
to the absorbing top electrode was reflected by a 
reflective top electrode. It seems to be resonant cavity 
effect caused by reflective top electrode. Also, it was 
known that light extraction efficiency of LED was 
increased by resonant cavity effect caused 

Fig. 3b shows the angle of incidence (0) and 
reflectivity (R) of an AlAs/AlGaAs DBR, which is 
used as the reflector for the reflective top electrode, for 
the different periods. In Fig. 3b, ~ 70% R and 15° 0 
were observed at the DBR of 4 periods. The 
reflectivity and angle of incidence were increased or 
decreased respectively, with increasing periods of the 
DBR because they have trade-off relationship [4]. At 
12 periods, ~ 70% R and 7.5° 8 were observed. 

For further 
light-current-voltage (L-I-V) characteristics of 590 nm 
AlGaInP LED chips with different structures were 
measured under an injection current ranging from 0 
mA to 40 mA as shown in Fig.4. In the current-voltage 
(I-V) characteristics, all of the samples show nearly the 


improved performance, the 


same operating voltages. The operating voltage of the 
LED chips developed was barely influenced by the 
presence and thickness of the TDBR because the AlAs 
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Fig. 4 (a) Comparison of the L-I-V characteristic curves of 


the conventional AlGaInP LED chip and AlGaInP LED chip 
with RTE having the TDBR of various periods. 


and AlGaAs materials used for the reflector have much 
lower resistivity than the underlying structures [13]. 

Distinct results were observed in the L-I 
characteristics of the 590 nm AlGaInP LED chips with 
different structures. Clearly the LED chips with a 
reflective top electrode show a higher output power 
than the conventional LED chips. Also, it was found 
that the higher output power was obtained at DBR of 
12 periods. At a high injection current of 30 mA, 
output power of | mW obtained for a DBR of 12 
periods shows relative increase of 23% and 18% 
respectively, as compared to without and with 4 
periods. The remarkable difference of output power 
between LEDs without and with 4 periods was not 
observed. Therefore, this results presented that cavity 
effect caused by insertion of RTE was depends 
strongly on reflectivity of DBR used. 


4. Conclusions 


In this work, efficiency improvements for a 590 nm 
AlGaInP light-emitting diode (LED) with a reflective 
top electrode (RTE) were investigated. A distributed 
Bragg reflector (DBR), consisting of AlAs/AlGaAs 
pairs, was used as a reflector for a reflective top 
electrode. It was found that a higher output power was 
observed from the AlGaInP LED with a RTE as 
compared to a conventional one. The output power of 
the LEDs with a RTE depends strongly on the 
reflectivity of the DBR used as its reflector. At an 
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injection current of 30 mA, the highest output power of 
1 mW was obtained from a 590 nm AlGaInP LED 
having the 12 pairs DBR with ~85% R, which showed 
a relative increase of around 23% over a conventional 
LED. 
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Abstract: The compensation method aiming to reduce the irradiation on the user’s body and to form a weak field area (a dark spot) in 


an antenna near region is considered. The results of the method development for the use in a narrow frequencies range are given. They 


show that under realistic conditions the compensation method allows reducing the irradiation power by a factor of 3-4 and more. The 


pattern changes are predictable and small, if the points of antennas placement are selected correctly. The dark spot necessity near the 


transmitter is not limited by mobile communication. A man often uses a low-power radio station, which is placed nearby in a vehicle, 


and it is necessary to protect from irradiation not only user but also other passengers. In the paper, it is shown that broadband field 


compensation can be achieved in various structures where the needed anti-phase fields are created either by identical antennas or by 


radiators located at equal distance from the compensation point, or by flat reflectors. 


Key words: Antenna theory, compensation, mobile antennas, mutual coupling, near fields. 


1. Introduction 


The compensation method was proposed to protect 
user’s head against strong electromagnetic fields of a 
cellular phone antenna. The method was intended to 
decide the particular problem of decreasing near region 
field created with the antenna of the certain device 
operated at the one given frequency. A cellular phone is, 
in particular, such a device. During a phone 
conversation, it is placed next to user’s head and 
irradiates sensitive human organs (brain, eye, etc.). The 
absorbed energy in today’s cellular phones can be as 
high as a half of all radiated energy. 

But the problem is not limited to the cellular phones, 
since, firstly, a man often uses or services a low-power 
transmitter, which he carries or places nearby in a 
production area or in a vehicle. Secondly, one must 
protect against a strong electromagnetic radiation not 
only living organisms, but also devices. The protection 
of devices is necessary, since the irradiation of nearby 


devices can disrupt their normal operation, cause 
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spontaneous device switching on and switching off, 
change operating regimes, etc.. 

The rather obvious idea of head protection by means 
of screening, i.e., by shading effect, is not practicable. 
The near field has no ray structure, and hence the 
shadow formed behind a metal screen can only cover 
an area approximately equal to the screen size. For 
example, in order to protect the head of the cellular 
phone user, the screen must be much larger than the 
cross section of the handset housing. For similar 
reasons, one should discard the idea of using an 
absorber, i.e., a dielectric shield that absorbs some part 
of energy. The distortion of the antenna pattern is 
another obvious disadvantage of using screens and 
absorbers. 

The protective action of the compensation method 
proposed by Bank [1] is based on another principle: the 
mutual suppression of fields created by various 
radiating elements in a certain area. The diverse 
variants of the principle implementation (interference 
of two linear radiators fields, folded antenna, and 
cavity antenna [2, 3]) possess a potentiality to reduce 
irradiation of user’s organism, especially his head, 


120 Compensation Method: New Perspectives 


without distorting the antenna pattern in the horizontal 
plane. 

In accordance with the key variant of the 
compensation method, as shown in Fig. 1, main 
radiator 1 is supplemented with second (auxiliary) 
radiator 2 situated in the plane passing through the head 
center and the feed point of the main radiator. The 
second radiator is placed between the head and the 
main radiator and is excited approximately in 
anti-phase to it (not exactly in anti-phase because the 
phase of the field progresses along the interval between 
the radiators). So, the radiator fields will compensate 
each other at a certain point inside the head, and the 
point will be surrounded with a zone of a weak field: a 
dark spot will arise. 

In order to obtain the same field strengths at the 
compensation point, when the radiators driving 
currents differ substantially, it is enough to separate the 
radiators by a small fraction of the wavelength, e.g. a 
distance of about 1-2 cm at 1 GHz. The relatively small 
driving current of the auxiliary radiator ensures that 
upon appropriate placement of the radiators the total 
field at far region is almost unaffected by the presence 
of the auxiliary radiator, and the pattern remains 
omnidirectional, allowing to decide the crucial 
question in the development of any method of reducing 
exposure, how to preserve the communication 
performances of the system. 

Results of key variant analysis are described in 
Section 2. 


ZT 2 


Fig. 1 Two radiators next to the head. 


2. Compensation at One Frequency 


The problem of field calculation in the space 
surrounding a radiator is complicated due to the 
heterogeneous nature of this medium. Cellular handset 
antennas are in close proximity to the user’s head, hand 
and body, which contain multiple tissues with different 
permittivity, much higher than that of the free space. If 
an antenna is tangent to the user’s head and is 
considered in a crude approximation as a linear radiator 
located along the (plane) boundary between two 
half-spaces: air with ¢. = 1 and a head with ¢, # 1, then 
an analogy is valid for a quasi-stationary electrical field 
(including an antenna field in the near region) created 
by alternating linear current and an electrostatic field of 
charged conducting bodies. On the other hand, it is 
known that an electrostatic field of an isolated wire 
located at the boundary of two media coincides with its 
field in a homogeneous medium with a permittivity 
equal to the average of both media ones [4]. 

This analogy allows reducing the calculation of an 
antenna near field in a piecewise-homogenous medium 
to approximate field calculation in a homogenous 
medium. But the near field of an antenna has a 
quasi-stationary character only in the first order 
approximation. The radiators are located near the head, 
but not tangent to a head. An assumption that radiators 
are located along on interface of two half-spaces is too 
crude. Therefore, the results are in need of check-up 
with the help of calculation program, for example CST. 

The circuit of radiators displacement is given in 


Fig. 2. Here, A; and Ay are the main andthe auxiliary 


Fig. 2 The diagram of phone radiator placement close to a 
human head. 
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radiator, C is the compensation point, O is the head 
center. The radiators are the symmetrical vertical 
length. The field 
components of the main radiator 4; in cylindrical 


electric dipoles with finite 


coordinates (p, g, z) are given by: 


SOY 1_Jexp(— jkR,)/R, + 


EE. = . 
é, sin kL, 


+exp(— jkR, )/R, —2cos kL, exp(— jkR,)/R, | 
y= ig lle) Fak) + 
(1) 

+(2+L,)exp(—jAR,)/R, —2zc0skL, exp(—jkRy)/ Ry | 
here, J4; is the driving current, Z; is the dipole arm 
length, R, =./(z—L,)’ +p; is the distance from the 
observation point M to the radiator’s top edge, 
R,=(z+L,)' +; is the distance from M to the 
radiator’s bottom edge, R, =|? +p? is the distance 
from M to the radiator center, and p, =p. For the 
auxiliary radiator A, one can write analogous 
expressions with L = Lp, pr = pr + by - 2p\bcosppr’, 
where, b is the distance between antennas. The field 
components of a monopole of the same arm length are 
smaller by factor 2. 

At the compensation point C, the both radiators 
fields must have the same magnitudes and opposite 
signs. But (1) shows that, in general, simultaneous 
cancellation of £, and E, by adjustment of the current 
J42 is impossible. Since, EZ, is greater than £,, and EF, = 
0 along the ¢ -axis (at z = 0), we choose to cancel the 
E, components: 

E.,(b+ ,,0,0)=—E.,(b + p,,0,0) (2) 
here, o is the distance from antenna A» to the 
compensation point. 

For the sake of simplicity, we assume that the 
dipoles have equal dimensions. The case of dipoles 


with different lengths is studied in Ref. [5]. 
The mutual coupling between closely located 
results in variations in their 


antennas some 


characteristics. Find the current and the input 


impedance of two adjacent finite length monopoles, 
located within the near region of each other [6]. In 
accordance with the theory of folded radiator and the 
superposition principle, one can analyze the 
two-radiator system with the generator exciting one of 
the wires as a superposition of two sub-systems: an 
even-mode sub-system (with in-phase currents), and an 
odd-mode sub-system (with anti-phased currents). The 
terminal currents Jy; and J49 of the radiators excited by 
the generators are e€, and @, are 

J 4 =(e, —e, JY, +(e, + é, JY, 

J 42 =(e, —e¥, +(e, +e, )Y, (3) 
where, Y; = 1/[(-;WicotkL + 2R)] is the input 
admittance of the long line formed with two monopoles, 
Y) = 1/((4 Z,,(L, a.) + 2R)] is the input admittance 
depending on the admittance Z,,(L, a.) of the monopole 
of the length Z and radius a, = Jab, W, = 120In(b/a) is 
the wave impedances of the line, & = 2n/A is the 
propagation constant, 7 is the wave length, R is the 
output impedance of each generator, b is the distance 
between the monopoles, and 2a is its diameter. 

Accordingly the self Z,, and the mutual Z, 
impedances of the radiators are: 
Zig 25 = 0.25 (1/ 1/7, 

=Z,,(L,a,)— j0.25W, cot kL + R, 

Z,) = 0.25(1/Y, =1/¥,) 
=Z,,(L,a,)+ j0.25W, cot kL.. 


(4) 


It is easy to show that in order to nullify the total 
field at the compensation point ( p, + 5,0,0 ), ie., 
E, (0) = E,(0)+ E,,(0) =0, the ratio of the radiators’ 


emf’s must be equal to: 
@ _ (+h )Fot (NF 6) 
&  @-YH)Rot¥ +h )ho 


here, F',9 is the value of function /F,, at the 


compensation point, where, F,,, is dependence of the 
field of the current excited with each generator from 
the distance to the observation point. Eq. (5) allows 
finding the voltage amplitude and phase at the input of 
the second radiator, if the amplitude and phase at the 
input of the first radiator are known: 
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J plI a = —Fig/ Fay 
i.e., the expression for required current ratio coincides 
with this expression in the absence of mutual coupling 
between the radiators. This fact occurs when the 
radiators are of the same length. For Z = A/4, we obtain, 


in the first-order approximation: 


J 9 Ry ° 
cm Re e pl J ( 10 20 )b 


where, Ry = Ip? + pe, Ry = 


As it follows from the Eq. (6), the dipole moments of 


LD? +(py +b)’ 


linear radiators must be in inverse proportion to the 
distances from the radiators to the compensation point. 

Fig. 3 shows two variants of the dark spot position 
relative to the head, where the spot boundaries are 
shown by a dotted line and the solid curve shows the 
head boundary. The feed points 4; and A> of the main 
and the auxiliary radiators and the compensation point 
C in Fig. 3 lie on the structure axis, and hence the dark 
spot is symmetric relative to this axis. In Fig. 3a, the 
compensation point is placed outside the head, and in 
Fig. 3b it is inside the head, close to its surface. In the 
latter case the whole dark spot lies inside the head, 
resulting in a significant reduction of losses inside the 
head. 

The dark spot dimensions can be calculated using 
equations similar to Eq. (1) written for the radiator’s 
field. The dark spot boundary is determined by radius 
Pn depending on b, z and g, at which n = |E/Eo| is 
equal to a given value 7 (e€.g., 9 = 0.1), where E, is the 
total field and £,,9 is the field of the main radiator. The 
intersection points , and > of curve n with a given 
level no determine the dark spot boundaries in the p 
direction and the length p;) of the segment, along 
which n is smaller than the required value mo. The 
length Ap= p2— p, of this segment, i.e., the dark spot 
length, is presented in Ref. [6] in the table and figures 
for different values of b, g and no, i.e., for different 
levels of field decrease at the dark spot boundary. 

The calculations show that upon optimal design of 


the structure, in terms of the relative positions of the 
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Fig.3 The dark spot: compensation point is located outside 
(a) and inside (b) the head. 


radiators and the compensation point, the volume of the 
dark spot increases, and the field inside this spot 
substantially decreases. One can choose the structure 
parameters so that the spot dimensions may coincide 
with the human head dimensions, resulting in a big 
reduction of losses power in the head, and in a high 
tolerance with respect to the user’s head movements. 
In Ref. [7], it is shown that using a second vertical 
auxiliary radiator placed along the symmetry axis of 
the structure, i.e., on the horizontal straight line passing 
through the head center O and the compensation point, 
permits to increase the dark spot volume. Employing a 
few of auxiliary radiators allows a further increase of 
the dark spot volume. At the same time the irradiation 
reduction factor increases. This factor is defined as the 
ratio of the losses power in the absence and presence of 
the auxiliary radiator within the limits of the dark spot. 
It is assumed, for the sake of simplicity, that the 
electric field strength increases linearly from the 
compensation point ¢ = 0 to the dark spot boundary t = 
to such that | £./E.19| = ns, where s = t/t, n=const(s), 
and hence the losses power exhibits a quadratic growth 
with the distance. The total losses power inside the spot 


is given by: 
1 
P =P, (ns)'ds =Pyn’/3 
0 


where, Po is the losses power due to the main radiator 
in the absence of the auxiliary radiators, within the 
limits of spot. For mo = 0.2 the losses power is 
diminished by a factor of 3/(0.2)° = 75 (or 18.8 dB). For 
No = 0.1 the reduction factor is 300 (24.8 dB), and for 1 
= 0.04 the factor is 1,875 (32.7 dB). In practice, this 
factor is roughly half of the calculated value, since the 
field within the dark spot has a more complicated 
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structure, i.e. the reduction factor is approximately 
equal to 1.5/no°. 

The pattern of the structure is considered in Ref. [8]. 
As said before, in order to preserve the circular form of 
the pattern, and in particular to avoid a deep gap along 
the structure symmetry axis, the auxiliary antenna field 
must be small compared to the main antenna field, i.e. 
this antenna must be placed at a minimal distance from 
the compensation point. If J; is the driving current of 
the main radiator, E., is its field, and J4, = Jy.De!” is 
the driving current of the additional radiator, where D 
is the ratio of the radiators’ dipole moments, then the 
field of the additional radiator is E)= E,De, where, u 
= kl — y,1is the path-length of rays from the radiators 
to a point located at angle @, and yw is the phase 
difference between the currents of the main and 
additional radiators. 

Since the total field of both radiators is equal to E = 
E, + E>, one can write for the pattern in the horizontal 
plane: 


|F| =|1+ Dexpl ju)/[(1+ Deosu, )° + D sin? uw, (7) 


here, uw, is the value of u at the denominator maximum. 
Since, D and u,, are constant values, the denominator of 
this expression does not depend on 4Q, i.e., the ratio 
between the pattern maximum to the minimum is equal 


to: 


_1£D (8) 


IF, 


min 


where the upper sign corresponds to a positive D and 
the bottom sign, to negative. It is easy to show that 
when the pattern maximum and minimum differ by 6 
dB (deviation of 3 dB from the average level), we 
obtain |D| = 0.33, and when the difference between 
the maximum and the minimum is 3 dB, we get |D| = 


0.17. 

Taking in the first approximation that Ri = 9 + b, Ro 
= £o, we obtain D = Ro/Ri = Po/(Po + 6), which reduces 
to Po = b/f, where f= 1/D — 1. Consider the case when 


the main radiator is situated on the structure axis, at a 
distance of 3 cm from the head, i.e., 9 + b = 3, and the 
compensation point is placed on the head surface. In 
this case, if D = 0.33, then # = 2, po = b/2 = 1, b= 2; 
for D = 0.17 we have £ = 4.88, 9 = 0.5, b = 2.5. The 
results point out the necessity for a tradeoff between 


the irradiation reduction level and the radiation pattern 
of the two-antenna structure, especially in cases when 
the antenna platform is subject to spatial restrictions. 
However, the pattern changes are predictable and small, 
if the points of antennas displacement are selected 
correctly. 

The simulation performed with the CST program 
show that the results of calculations, simulations and 
measurements are in good agreement, if the user's head 
is absent. The coincidence between the results in the 
presence of head is not too good (in accordance with 
the approximate estimation of permittivity), but a 
qualitative agreement does exist. 

The obtained results show that under realistic 
conditions the compensation method allows reducing 
the irradiation power by a factor of 3-4 and more. 

A method of user’s irradiation reduction similar to 
the compensation method is proposed in [9], where the 
concrete ratio of two monopoles driving currents 
required for decreasing the SAR level at 2 GHz is 
found by calculations using the Method of Moments. 
But in this paper an analysis of causes of an irradiation 
reduction and recommendations for the use of the 


observed effect are absent. 
3. Wideband Field Compensation 


Cellular communication is not the only area of 
application, where it is necessary to create a region of 
weak field near the transmitting antenna. The task is 
vital one, if, for example, a man uses or services a 
low-power radio station, which he carries or places 
nearby in a production area or in a vehicle, where it is 
necessary to protect from irradiation not only user but 
also other workers or passengers. Creating a weak field 
region in such cases is an efficient technique of 
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protecting against irradiation. 

In these cases, requirements to the produced weak 
field change in part. On the one hand, they often get 
complicated due to necessity of operation in a broad 
frequency range. On the other hand, the requirements 
become less stringent. For example, large dimensions 
of the transmitter allow increasing the room for 
mounting auxiliary radiators. Presence of a signal in 
the direction opposite to user’s location sometimes is 
not required. The requirement for the horizontal pattern 
to be flat in the direction of user’s location is somewhat 
lifted. 

Operation in a broad frequency band gives rise to 
additional difficulties, since the fields of the main and 
the auxiliary radiators at the compensation point must 
have the same magnitude and opposite sign at each 
frequency. The phase shift can be easily implemented 
by using a transmitter output stage with balanced 
output. But, firstly, these two instances of the radiator 
type need to have identical characteristics at all 
frequencies of the operation range. Secondly, since the 
main and the auxiliary radiators are located at different 
distances from the compensation point, and signal 
velocities in the air and in the feed lines connecting the 
radiators to the transmitter are different, it is necessary 
to insert a small phase shifter in one channel. This 
phase shifter must provide the required delay time 
throughout the frequency range. 

It is not easy to meet both these requirements. To 
provide the same signal amplitudes, it is necessary to 
secure identity of the main and the auxiliary radiators 
or to insert controlled attenuator in one channel. 
Wideband antennas often have a complicated structure, 
which uses different conducting and _ dielectric 
materials. That is why antenna characteristics are 
unstable, i.e., different antenna instances have different 
frequency responses, and so their fields are not equal. It 
is impossible to construct controllable attenuators in a 
wide frequency range. 

Placing both radiators at equal distance from the 
compensation point yields better results. 


A simple structure consisting of two radiators placed 
at the same distance from the compensation point is 
shown in Fig. 4. The radiators are two asymmetrical 
antenna instances of the same type, which have the 
circular pattern in the horizontal plane and similar 
amplitude-frequency responses. They are located at the 
same distance b = a/cosf from point C, where a = AC 
is the distance from point C to the middle line of the 
transmitter cover. 

For example, the auxiliary radiator phase differ from 
the main radiator phase by 180°, and the amplitude of 
the vertically polarized signal £, of the auxiliary 
radiator be £, = DE,, where £, is the amplitude of the 
main radiator signal, and D < 1. If D = 1, ie., the 
amplitudes of the auxiliary signal and the main signal 
coincide, the sum of the signals in the plane of the 
structure symmetry is zero, since the signals in the 
plane are the same in amplitude and opposite in phase. 
If D # 1, the summing signal EF = £, + Fy = (1 — D)E; 
is equal to a fixed small share of the main signal and 
weakly frequency dependent, since the distance from 
the radiators to the compensation point is the same, and 
the radiators themselves are similar. No phase shifter is 
needed in this case. 


Fig. 4 The two-radiator structure. 1 is the transmitter top 
cover, 2 is the user’s body, B, is the main radiator, B, is the 
auxiliary one, and C is the compensation point. 
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A disadvantage of the proposed circuit is the absence 
of the far field signal in the structure symmetry plane 
(in both directions), when £» = £;. But, first, this is not 
critical in many applications. Second, the angle (gap 
width), in which radiation signal is nearly zero, is very 
small. And third, the influence of the ground and 
neighboring bodies results in a smoothing of the pattern, 
i.e., both the depth and the width of the gap in the 
horizontal pattern will be reduced in real conditions. 

Estimate the gap width by calculating angle Ag 
between its boundaries, i.e., between the points 
corresponding to given small value /; of the pattern. If 
a = 25cm, and f = 30°, the signal at angle g, which 
consists of two signals with equal amplitude /, and 
opposite phase created by radiators located at distance 
d = 2atanf = 28.9 cm from each other is equal to E = 
2E\sin[(kd/2)sing], i., Emax = 2E,. Here, the 
cylindrical system of coordinates is used (Fig. 4). If @ 
is small, 

f, =sin|(kd/2)sin gp] ~ kde/2 (9) 
Le., Ag = 29 = 4f,/(ka). In particular, if f; = 0.1, the 
wavelength 2 is 30 cm, and d = 28.9 cm, then Ag = 
3.8°. 

Another disadvantage is the fact that, when D # 1, 
the irradiation power reduction factor inside the dark 
spot decreases, as the total field, fails to vanish even at 
the compensation point. But the disadvantage is also 
inherent in the case of antenna placement at different 
distance from the compensation point, due to the 
limited dimensions of the transmitter top cover. 

Experimental results for the signal magnitude at the 
compensation point are given in Fig. 5. Curve 1 
corresponds to the signal of the main antenna placed at 
point A (Fig. 4), and curve 2 corresponds to the signal 
of two antennas placed at points B, and Bo. As seen 
from the results, mounting both radiators at the same 
distance from the compensation point permits to ensure 
compensation in a wider frequency band. The auxiliary 


radiator decreases the signal magnitude at the 


compensation point by 10-15 dB (3-6 times) in the 
band from 1.7 to 2.7 GHz. 

Using of two auxiliary radiators gives additional 
advantages. In this case, the amplitude of each 
auxiliary radiator signal may be smaller than that of the 
main radiator. One variant of this compensation 
structure is presented in Fig. 6 (top view) where the 
same notation is used. 

If Eo is the amplitude of the main radiator signal at 
the compensation point, the amplitude of each auxiliary 
radiator signal must be EF, = £y = E)/2. In the 
cylindrical coordinate system, the origin of which 
coincides with point A, the ratio of the auxiliary signal 
to the main signal in direction gis: 

E,(p)/E,(e)= 0.5 exp Lik (a, —b,)h 

E,(p)/E,(e)=0.5exp[jk (a, +5,)] (10) 
where, Fo() is the main radiator field at distance p, 
E;(p) is the field of the auxiliary radiator located at 
point B;, and a, + b, are the path-length differences of 
the signals from the main and auxiliary radiators to the 
observation point. Here, a; = a(1 — cosf)sing, b, = 
asinfcos¢@ (these segments are marked in Fig. 6). The 


total signal is: 
E(p)= E,(p)- cos(kb, exp(jka,)} (1) 
The ratio of the total and the main signals is equal to 
E(p)/|E.(p) = [1 - cos (kb, Joos (ka, JF + 


+cos?(kb, )sin* (ka, ) (12) 


In the structure symmetry plane, the total signal in 


the far region is: 
E(p)/E,(p) =1-exp| jka(1—cos £) | = 
= -2j exp| jka(1—cos B)/2 |sin| ka(1—cos B)/2 | 
Le., 
|E(e)/E,() =2sin[ka(l—cosf)/2] (13) 

If, for example, a = 25 cm and £ = 30° or 60°, we 
obtain at 1 GHz |E(py/Eo(pl = 0.69 and 1.22, 
respectively. The total far field in the symmetry plane 
in both directions is other than zero. It follows from 
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Fig. 5 The frequency dependence of the signal at the 
compensation point, when the main radiator (1) and two 


radiators (2) are placed at the same distances from the point. 


Fig.6 The three-radiator structure. 
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Fig. 7 The horizontal patterns of different structures at 1 
GHz. 


Eq. (6) that at any g the total signal in the far region is 
other than zero, too. This is an inherent advantage of 
such structure. 

The calculated horizontal patterns at frequency 
1 GHz of one radiator placed at point A, of two 
radiators (Fig. 4) and of three radiators (Fig. 6) located 


on the cover of infinite dimensions are given in Fig. 7 
and marked with symbols 1, 2 and 3, respectively. 

The considered compensation structures comprise 
radiators placed at the same distances from the 
compensation point. The structures require no use of 
any phase shifter. The compensation structures, in 
which flat metal reflectors (mirrors) are used instead of 
auxiliary radiators, permit to create a weak field in the 
subspace on the side of the reflectors, where the 
antenna is located, which is similar to the one generated 
by three radiators. 

An example of such structure is presented in Fig. 8 
(top view). As seen from the figure, the flat metal plate, 
in which main radiator A is reflected, is mounted on the 
transmitter top cover along the structure axis of 
symmetry (in the vertical plane passing through 
compensation point C). The reflection (or imaging) of 
the main radiator signal by the metal plate into the 
space in front of the reflector is equivalent to the 
presence of second radiator B located at the same 
distance behind the plate. 

The reflected signal phase differs from the main 
radiator phase by 180°, since vector £, of the incident 
signal is located along the tangent to the reflector 


surface. 


Fig. 8 The structure with one reflector. R is the metal 
reflector, and B is the location of the equivalent radiator. 
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The amplitude of the reflected signal E yis Ey = DE, 
where D depends on magnitude @ of the angle, 
subtended at point A (D < 1) by the reflector. 

If the plate dimensions are infinite (or if they are 
much larger in all directions than the main radiator 
dimensions), then D = 1, i.e., the amplitudes of the 
reflected and the main signal coincide, and the total 
signal in the reflector plane is zero. If D # 1, total 
signal F = £, + Ey = (1 — D)E;, isa fixed share of the 
main signal, and the share is independent of frequency, 
since the distances from both radiators to the 
compensation point are the same, and the radiators 
themselves are identical. 

On the whole, the performance of the structure is 
similar to that of the two-radiator structure presented in 
Fig. 4. No phase shifter is needed in this structure. An 
additional and very important advantage of the 
described structure is that the external actions, e.g., an 
approach of metal bodies, exert practically no effect on 
the structure operation, as the main radiator and the 
reflected signals undergo the same changes upon any 
external actions. Only the appearance of a metal body 
between a signal source and the compensation point is 
an exception and a very rare one at that. 

The influence of a mirror (a metal plate) on the 
antenna input impedance decreases the matching 
between the antenna and the transmitter. To offset the 
decreased matching, it is necessary to increase the 
distance between the antenna and the mirror as well as 
the mirror width so that the second signal would not 
decrease. The reflector width has only a weak effect on 
the antenna input impedance, but essentially increases 
signal E,. Moreover, the mirror influence on the 
antenna input impedance is constant and can be 
allowed for a priori. The mirror height must exceed that 
of the antenna. 

Flat reflectors can be manufactured in the form of a 
light-weight, strong and collapsible construction (if 
necessary). In such a case, each reflector is 


implemented as a set of vertical parallel wires, located 


at a distance of 0.06A nin, where, Amin is the minimum 
wavelength in the frequency range. 

Using two reflectors, i.e., creating two equivalent 
radiators, gives additional advantages. In this case, the 
amplitude of each reflected signal may be smaller than 
that of the main radiator signal, and so the requirements 
to the reflector loosen. A variant of such compensation 
structure is presented in Fig. 9 (top view), where, as 
before, B, and Bj denote the equivalent radiators 
locations. The performance of the structure of Fig. 9 is 
similar to that of the three-radiator structure presented 
in Fig. 6: the distance from all radiators to the 
compensation point is the same, i.e., field variations 
due to frequency dependence have no effect on the 
compensation quality. 

As stated above, the amplitude of the reflected signal 
(vertically polarized) is Ey = DE). Here, D depends on 
the mirror dimensions. 

An essential drawback of the structure with two 
reflectors is that the reflectors creating the dark spot in 
the near region block at the same time prevent the 
electromagnetic waves in the angular sector, shadowed 
by them. Here, the signal in the open sector increases 


Fig.9 Compensation structure with two reflectors. 
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This drawback limits the 
application of such structures. 


approximately twice. 


The structure with two reflectors was tested 
experimentally. A general view of the experimental 
setup is presented in Fig. 10. The following notation is 
used in the figure: 1 is the radiator, 2 is the receiving 
antenna, placed at the compensation point, and 3 and 4 
are flat reflectors, located symmetrically on either side 
of the straight line passing through the compensation 
point and the radiator. The distance between the 
radiator and the receiving antenna is 25 cm. The 
reflector is manufactured in the shape of a vertical 
rectangular plate. Contact with the ground is 
maintained along the whole width of the reflector, 
using a horizontal element. In the course of the 
measurements, each reflector was separately moved 
along a straight line to the point where the received 
signal at 2 GHz (close to the band center) decreased by 
6 dB. As a result the antenna field at the compensation 


point in the presence of both reflectors is substantially 


Compensation Method: New Perspectives 


smaller than the field of the antenna without reflectors. 
The measurements are performed at two different 
values of reflectors width b, i.e., 8 and 15 cm, and at 
different angles (10°, 15° and 20°) between the axis of 
the setup symmetry and the lines, along which 
the reflectors were moved to attain the zero field at the 


Fig. 10 General experimental setup. 


Table 1 The measurement results for the signal magnitudes at the compensation point, in dB. 


With reflectors of width (cm) Difference 
f (GHz) Without reflector 8 15 8 15 
B/2=10° 15 15 10 15 15 
0.5 -17 -32.1 -33.3 -35.8 15.1 16.3 18.8 
1.0 -19.5 -33.5 -35.8 -34.3 14 16.3 14.8 
1.5 -39 -46.4 -43.4 -54.9 7.4 44 15.9 
2.0 -28.8 -37.3 -40.4 -62.5 8.5 11.6 33.7 
2.5 -27.7 -47.1 -42.6 -41.3 19.4 14.9 13.6 
3.0 -36.3 -43.9 -47 -42.6 7.6 10.7 6.3 
a) E,aB b) £,dB 
-10 -10 
P yon 
-30 30 jaa 
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Fig. 11 Frequency dependence of signal at the compensation point without reflectors (1) and with reflectors (2) for reflector 


width (a) 8 cm and (b) 15 cm. 
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compensation point. In all variants, when the distance 
between the compensation point and the near edge of 
the reflector was 9 cm, the signal was half as strong as 
its value without the reflector. 

The measurement results are given in Table 1. The 
frequency dependence of the signal at the 
compensation point without reflectors (1) and with 
reflectors (2) is presented in Fig. 11, for reflector 
widths of 8 (a) and 15 cm (b). As it is clear from the 
table and figures, the reflectors with width 8 cm 
permit one to decrease the signal magnitude without 
tuning by 10-15 dB for the whole frequency range, 
and the reflectors with width 15 cm perform even 
better. 

The dark spot boundaries of the three radiator 
structure (Fig. 6) in the horizontal plane (z = 0), which 
are determined by the procedure described above, show 
that these spots have sufficiently large dimensions. 


4. Conclusions 


We have analyzed a compensation method for 
wideband operation. Different structures for wideband 
compensation are introduced and investigated, for 
example, replacing the phase shifter by a delay line. 
The calculation and experimental results show that using 
several radiators located at equal distances from the 
compensation point or using flat reflectors, one can 
decrease the irradiation power by a factor of 5 or more in 
a large bandwidth. Two auxiliary radiators or the flat 


reflectors permit to achieve these results without 
introducing deep and wide gaps in the horizontal pattern. 
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Abstract: To solve the u problem of the Minimal Supersymmetric Standard Model (MSSM), a single field S is added to build the 
Next Minimal Supersymmetric Standard Model (NMSSM). Vacuum enlarged with non-zero vevs of the neutral-even CP is the 
combination of H,, Haand S. In the NMSSM, the Higgs sector is increased to 7 (compared with 5 hogs in the MSSM), including three 
Higgs-which are the even-CP hy» 3 (Myi< My2< Mp3), two Higgs-which are odd-CP a,» (m,, < m,2) and a couple of charged Higgs H". 
The decays Higgs into Higgs is one of the remarkable new points of the NMSSM. In this paper, we study the decays H~ into W* and 
a,. The decay width is calculated to one loop vertex correction. The numerical results are also described together with evaluations. 


Key words: Higgs, decays, NMSSM. 


1. Introduction 


The simplest version of supersymmetry is the 
Minimal Supersymmetric Standard Model (MSSM). 
This version is limited by two problems: the p and the 
hierarchy [1-4]. The simple supersymmetry, which is 
beyond the MSSM, is_ the 
Supersymmetric Standard Model (NMSSM). The 
special characteristic of Higgs boson in the NMSSM is 


Next Minimal 


the decay of Higgs into Higgs. It is remarkable that the 
lightest state al of the odd-CP Higgs (CP is Charged 
and Parity transition) can play a role of a pseudo, 
goldstone, which has small mass and can lead to the 
predominated decay of the even-CP h— aja; [5]. The 
even-CP Higgs and the heavy odd-CP Higgs can be 
generated at Large Electron—Positron Collider (LEP) 
in e*e — ha, but they may not be discovered because 
the dominant h decay was not searched for. There are 


different ways to make the mass of Higgs boson 
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increased in the MSSM and in the beyond MSSM. One 
simple way is to study the beyond singlet of the MSSM 
which term AS H,H, in the 


super-potential, this is the term that contributes 


contains one 


A*v’ sin? 28 at v = 174 GeV to the squared mass of 
even-CP Higgs [6] and therefore, it can make the mass 
of Higgs boson increased over the limit of independent 
decay state. It should be noted that this contribution is 
maximum with tanB ~ 1. Thus, a condition in which the 
lightest odd-CP Higgs with its mass is under 2m,, the 
two lightest even-CP Higgs boson and the charged 
Higgs boson can be found in the MSSM. All of them 
can be generated at LEP and they are now being 
searched for. 

The charged Higgs makes up more than 40% in the 
top-quark decay at Tevatron; the products of this decay 
are charged Higgs and bottom-quark (t > H*b ). The 
decay method of charged Higgs is H* > W*a,, with 
a, > CC, gg. 

The neutral Higgs sector in the NMSSM includes the 
following states: three even-CP and two odd-CP. Many 
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analysis on Higgs sector in the NMSSM [7] have 
shown that, in the specific physical state of the even-CP 
Higgs, there is a strong mix between the doublet state 
and the singlet SU(2) with the reduction in the 
interaction of gauge boson. The study on light Higgs 
contributes to the discovery of one or more Higgs states 
at LEP, at Large Hadron Collider (LHC) [7] and at 
large energy accelerators. 

The NMSSM is established from the MSSM when 
the MSSM is added with a gauge chiral single 
superfield S, which includes the renormalized 
superfield interaction and the soft supersymmetry 
breaking term Ls. In the NMSSM, the terms of the 
super-potential Wyigg, are dependent on superfield 
Higgs H,, H, and S (here, we follow the SLHA2 
regulations, however A, is also written as A, and 


H,, is also written as H,): 


W, =(u+AS)H, -H, +E pS+uS? +2 9°) 


Higgs 


with: 4, « is the non-dimension coupling Yukawa; 
u, ’ is the supersymmetry mass; &, is the square 
supersymmetry mass parameter. 

From Eq. (1), Yuakawa interaction of quark and 
lepton superfield are added to: 

Wikia’ h, HH, -OUS + h,H OD; + hg LES (2) 
here, the Yukawa interaction h,, hg h. and the 
superfields O, US, De, 1. Es are the matrices and 
vectors in the corresponding spaces. 

The soft breaking supersymmetry sector is regulated 
in SLHA2: 


=a see 2: 
-L ‘soft ~~ Min, 


veh 


u 


2 2 2 2 
+m; |H,| +m; 


S|’ +m |0"| 


2 2 2 
mi |Uz| +mp 


D,|-+m;|L?|+ mz |Ep| : 
+(h,A,O.H,Up —hyAjO-H Dy -h, AL. ER 2 


1 
+AA,H,.H,S+ : KAS? +m3H,.H,+m’S’?+ESthe) 


(definition m; = Bu,m? = Bi’). 


In the super-potential Eq. (1), we _ have 
supersymmetry parameters 4, wv” and &, (the soft 


supersymmetry interactions break the parameters 


m;, mé and parameter &s in Eq. (3)), however, some 
terms are not eliminated in some different solutions for 
simple NMSSM with a part of invariant super-potential 
when p=p'=€, =0 


Woussy = ASH Ha +58 (4) 


Then, eliminating the parameters m;,m? and és in 
Eg. (3), combining the vevs of Sin the weak sector or 
in the breaking supersymmetry to define wu: 

Hage = AS (5) 

The matter of p. in MSSM has been solved then. 

As any supersymmetry theory with invariant 
super-potential sector (ternary), the Lagrangians, 
which contain the soft supersymmetry violation 
conditions specified by Eq. (3), have one Z3 symmetry 
randomly, which is corresponding to the multiplication 


2/3 The invariant 


of all chiral superfields with e 
super-potential Eq. (4) is presented like invariant Z3. 
The non-dimension terms in the super-potential Eq. (1) 
will break the symmetry Z3. The model with 
super-potential Eq. (1) is the NMSSM. The invariant 


Z3 Higgs sector is defined by the seven parameters 


2 
A, a 


A, K, my ,m m:,A,,A,. The expressions of Higgs 
mass matrix in the invariant Z3 of the NMSSM show 
that invariant Z3 is obtained when: 

m, =ms =&,=p=p'=€, =0 (6) 


From the supersymmetry gauge interaction and soft 


supersymmetry breaking conditions, we obtain Higgs 


potential: 
V riggs = |AH IH] — HOH? ) + eS? + 2u'S + E-| 
2 2 ol? +? 
+(mz, +|u+AS| (\H; +|H; ) 


+ (my, +|u+ asl (as) +|ajf 


+ BL EEE i snc la ap 


ey 


2 
§2 | py + py 0* 0 py 0 |? 21e/2 
ore Jn tHe +A He +m; (S| 


+(AA, (HH) -HOH)S + 5K AS" 


+m (HH; - HOH! )+meS* +&S +he 
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where gi and g> present gauge interaction U(1) and 
SU(2). 

The even-CP Higgs sector has three mixed states, 
which are the real parts of H,,, Ha and S. The largest mass 
state h, which is nearly the same with SM, interacts with 
electroweak Gauge boson, this state has squared mass 


M,. given by [7, 8]: 
M; =M? cos? 2B+A°v sin’? 28 + rad.corrs + A, (8) 


The diagonal elements of matrix in the singlet sector 
are given by (with the assumption that s > v,,v, ) 


M2 =ks(A, + 4ks) (9) 


The term A, in Eq. (8) is derived from the mixture 
of doublet and singlet 
od SAE 2p) (10) 
Mi, -M3, 


where, M; is given by M; and has no mixed term. 
And the mass of charged Higgs boson is defined by 


2 
M2=M2,,+Vv [#2] (11) 


The process of the decay of Higgs into Higgs also 
contains other interests; that is the decay of one 
charged Higgs H~ into W~ + a, (or W*+h ) 
with a, > gg,cqt't , this is the channel that was 
missed in the process of searching for the decays 
<t—H"b at Tevatron. In the MSSM, the mass of 
charged Higgs is related to the mass of the particle W 
and odd-CP Higgs. 

mi. = mM. +m, +(SUSY loop) (12) 

Supersymmetry (SUSY) does not contribute much in 
this relationship, since m, is minor, we can see 
that m,.*m,,... The dominant decay process of 
charged Higgs is H* > W‘a,, which is nearly 100%, 
except when the mass of charged Higgs almost reaches 
to the mass of Boson W’. In that case, the decay 
H* > Wa, is still 
B(H* — W*a,) => 70% 


process dominant, 


The searching for charged Higgs which decays into 


W*a, is relatively easy at LHC, where quark top 
produces charged Higgs (z7—H*b ). Studying the 
decay of new particles in the Model will bring us the 
hope of finding out these particles as well as verifying 
the correctness of that Model [9]. In this paper, we 
study the decay H* — W“a, and calculate the decay 
width to one loop vertex correction. The numerical 
results are also shown in graph to evaluate the 
influence of m, to the decay width and to the lifetime 
of H*. 


2. The Decay H* > Wa, 


The amplitude which is calculated to tree level has 
the following form: 


M, =&, [i.2)-£(sin pus, +cos BUs\(a+4,)" . (13) 


The decay width which is calculated to tree level is: 


2 
M,| 4 
I, (H* > W*a,) = [Mo[ 4 af (14) 
16m), 
with 
: ae We P py? 2 2 9 
| = —“L (sin BU,, +cos BUS) (2m;, +2m;, —my) 
2 B20 \? 
+5 (sin BUI, + cos BUS, ‘; (oe =m , 
4 My 
and A = mi +m, 4 m;, 2m;,m;, 2mm, 2mm, 
a1 
Sa 
ay 


Fig. 1 
decay H* —> W*“a,. (a) Tree level; (b), (c) and (d) one loop 


Feyman diagram for correction SUSY-QCD in 


vertex correction. 
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The decay width which is calculated to one loop 
vertex correction has the following result: 


Mal 4 


I(H* >W*a,)= (15) 


16m}, 


with |,| =M;-M,+2Re(M;-6M,+M; 5M, +M; -6M,) 


M; -6M, = ae ee (sing-Uf, +cosB-U; ) (sing-UZ +cos-U;. 
70 
(om 2m? —m},) +—=-(m? m) [esas 
Ww 
2 2 ade a 
| (mi, + 3m, ma) +i 7 au =) Hae) 
Ww 


with Co(L) = Col imi, my, mj, 15, ,mp, My 


and C71) = Cri m;,mj,,m;,,ms,,m>,,m;, ) are the 
Pasarino-Vellman functions. 
The vertex elements A, B, in which S; = hj; S. =h»; 
$3 = hg and Py =a,; Po =a», 
A=-igmy (O; cos B+U%, sin B) 


-8""_(ys sinB+U5, cos f3).cos 28 
2cos Oy, 


Wie 
+ i—(vU%, +v,U%,) sin 2p 


V2 
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g af 2g” 
POY 
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with: Co(3) = Col my, ,m;, , Mj, My M5, Mp, ) 
and C)2(3) = Cy2(m,,,m;, ,m;,,my,,Ms,,Mp,)- 

To study the influence of the mass m,; on the decay 
process H* — W*a,, we use two set of parameters [7, 
10] for programming numerical calculation. 

The 1st parameter set: 2 = 0.8; x = 200; 4=0.1; my 
= 100 GeV; tanf = 3; sina = —0.58; A; = 6; A; = 486. 
From the results obtained, we find that the influence of 
Mi on the decay H* > Wa, is relatively significant 
(Figs. 2 and 3). 


Specifically, the influence of m, on the decay 
width and on the lifetime of H* in the decay 
H* > Wa, is relatively significant. When ma runs 
from 4 to 5, it can influence 2.9% to the decay 


width and about 4% to the lifetime of H*. With the 1st 
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Fig. 2 The influence of m,; on the decay width of the decay 
H* > Wa, with 2 = 0.8; x = 200; k = 0.1; my = 100 GeV; 
tanf = 3; sina = —0.58; A, = 6; A, = 486. 
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Fig. 3. The influence of m,, on the lifetime of H* of the 
decay H* — W’a, with 2 = 0.8; x = 200; k = 0.1; my = 100 
GeV; tanf = 3; sina = —0.58; A, = 6; A, = 486. 


parameter set, we obtain the results of the decay width 
H* - Wa, at about 0.035-0.043 (1/s), and that of the 
lifetime of H* at about 23.5-29 (s). 

The 2nd parameter set: 1 = 0.8; x = 200; & = 0.1; 
m,, =100 GeV ; tanf = 10; sina = —0.726; A, = 7; A, = 
492. We obtain Figs. 4 and 5. 

From the graph results with the 2nd parameter set, it 
can be seen that the influence of m,; is also relatively 
significant, its contributions to the results are even 
more than that of the 1st parameter set. When m,; runs 
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Fig. 4 The influence of m,; on the decay width of the decay 
H* + W*a, with 2 = 0.8; x = 200; k = 0.1; my = 100 GeV; 
tanf = 10; sina = —0.726; A, = 6; A; = 486. 
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Fig. 5 The influence of m,; on the lifetime of H* of the 
decay H* — W’a, with 2 = 0.8; x = 200; k = 0.1; my = 100 
GeV; tanf = 10; sina = —0.726; A, = 6; A, = 486. 


from 4 to 5, it can influence 6.7% to the decay width 
and about 2.9% to the lifetime of H*. With this 


parameter set, we obtain the value of the decay width 
H* > W‘a, at about 0.015-0.018 (1/s), and that of the 


lifetime of H* at about 55-61 (s). 
3. Conclusions 


In the NMSSM, a single superfield is added with 
complex scalar field components, this causes the 


appearance of seven Higgs in the NMSSM. In the 


H* Decays into W* and a; in the Next Minimal Supersymmetric Standard Model 135 


search for Higgs boson basing on decayed products, we 
need to focus on some channels which have high 
branching ratios. In this paper, we study the dominant 
decay of charged Higgs: H* > W~a,, it can be nearly 
100%, except when the mass of charged Higgs almost 
reaches to the mass of W’. 

The above decay width is calculated to one loop 
vertex correction. The numerical results which assess 
the influences of ma; are also shown in this paper. In 
this decay, the influence of m, is relatively significant. 
We need to pay attention to this element when studying 
theories as well as paying attention to the decay 
experiments of charged Higgs. The lifetime of H* is 
relatively short, and this brings us the hope that we can 


find charged Higgs soon. 


Acknowledgments 


This research is supported by the National Foundation 
for Science and Technology Development (NAFOSTED) 
of Vietnam. Grant number: 103.03-2012.80. 
References 


{1] D. Radovan, F.G. John, Many light Higgs bosons in the 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


next-to-minimal supersymmetric model, Phys. Rev. D 79 
(5) (2009) 055014-055029. 

M.M. Almarashi, S. Moretti, Very light CP-odd Higgs 
bosons of the NMSSM at the LHC in 45-quark final states, 
Phys. Rev. D 84 (1) (2011) 015014-015021. 

H.E. Haber, G.L. Kane, The search for supersymmetry: 
Probing physics beyond the standard model, Phys. Rep. 
117 (2-4) (1985) 75-263. 

W. Bernrenther, M. Suzuki, The electric dipole moment of 
the electron, Rev. Mod. Phys. 63 (2) (1991) 313-340. 
M.M. Almarashi, S. Moretti, LHC signals of a heavy 
CP-even Higgs Boson in the NMSSM via decays into a Z 
and a light CP-odd Higgs state, Phys. Rev. D 85 (1) (2012) 
017701-017705. 

A. Bartl, H. Eberl, K. Hidaka, S. Kraml, W. Majerotto, W. 
Porod, et al., SUSY-QCD corrections to stop and sbottom 
decays into W and Z bosons, Phys. Lett. B 419 (1998) 
243-252. 

U. Ellwanger, Higgs bosons in the next-to-minimal 
supersymmetric standard model at the LHC, High Energy 
Physics -Phenomenology (2011). 

U. Ellwanger, C. Hugonie, A.M. Teixeira, The 
next-to-minimal supersymmetric standard model, Phys. 
ept. 496 (1-2) (2010) 1-77. 

N.C. Cuong, Squark decays into H° in the MSSM with 
complex parameters, Commu. Phys. 14 (4) (2004) 
231-237. 

R. Dermisek, Light charged Higgs in the NMSSM, 
hep-ph/1012.3487 (2010). 


v3) 


Journal of Physical Science and Application 4 (2) (2014) 136-138 


PUBLISHING 


The Statistical Distribution for g-Deformed Crystal 


Lattice Vibration for Generic Atomic String 


Nguyen Thi Ha Loan 


Department of Physics, Hanoi Pedagogical University No 2, Vinh Phuc 280000, Viet Nam 


Revised: August 30, 2013 / Accepted: October 30, 2013 / Published: February 15, 2014. 


Abstract: In this paper, a g-deformation of crystal lattice vibration for generic atomic string is suggested. The expression for energy 
spectrum depends on g-deformed parameter have been obtained. The statistical distribution of -deformed crystal lattice vibration is 


calculated and this result is very much different from there of q-deformed crystal lattice vibration. 


Key words: Quantum, crystal, lattice, vibration. 


1. Introduction 


Over the last decades, there has been a lot of interest 
in the study of quantum groups and algebra. This is 
dictated by the point of view about their applicability in 
physical models [1-4]. They are relevant for a variety 
of diverse problems, such as the study of solutions to 
the Yang-Baxter equation, rational conformal field 
theory, two dimensional field theory with fractional 
statistics. 

The deformed oscillators formalism plays a useful 
role in the study of the quantum optics, condensed 
matter physics, etc.. 

In this work, we construct a §-deformed crystal 
lattice vibration for a generic atomic string; this is a 
fundament for study of crystal lattice vibration, which 
would give a base for a new quantization procedure. 


2. &-Deformed Oscillator 


The Hamitonian of &-deformed oscillator, system is 
given by [5, 6]: 


ps 


H =—+=—mo’x’ (1) 
2m 2 
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We define the creation (annihilation) operator a” (a) 
of &-deformed oscillator as: 


pes 70 x= P) 
2h mao (2) 


Hence, the coordinate and momentum operators can 
be expressed in terms of creation and annihilation 


+ 
operators a,aas: 


x= b (a + a‘) 
2m (3) 
p= «|m( - a‘) 
2 
Using Eqs. (3) and (1), we obtain: 
H="2(a-a" +a" -a) (4) 
2 
The creation (annihilation) operator a-(a) and the 
number operator N_ satisfies the commutation 
relations [2]: 
a-a’—g-a’-a=1 
[g.a =| ga" | =0 
=8*, g?=1 (5) 
N, a| =-a 
a a’ | =a’ 


The Fock space basics can be introduced: 
_@y) (6) 
= ~——/0 
y= ro) 
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where, |°) is the ground state satisfying the conditions: 
a| 0) =0 


N|n)=n|n) (7) 


The energy spectrum of the &-deformed oscillator is 
defined by the equation: 


n) (8) 


where 


(9) 


3. The §-Deformed Crystal Lattice Vibration 
for a Generic Atomic String 


The Hamitonian of generic atomic string can be 


expressed as: 
if 1 1 
H= Yaz PR +5 Mora) (10) 


ux 18 genezalized coordinate and P, is 


where, 
corresponding momentum of generic atomic string. 
x is total of k, which is waved vector in the first 
Brillouin area. 

annihilation 


The creation and operators 


corresponding to wave vector k have the forms: 


= peel Vion, + eh) 


corre Ji (11) 
es 1 te 


The coordinate and momentum operators can be 
expressed in terms of creation and annihilation 


operators @> % as: 


h sa 
u,=, Sage tee) (12) 


Using Eqs. (12) and (10), we obtain: 
1 
H=> "Pla a +aja,) (13) 


k 


The oscillator number operator and the oscillator 
creation, annihilation operator satisfies the following 


commutation relations: 


G30) oie ga" a, = Oey 
a,a’, — ga*,.a; = 0 (14) 
[g.a,]= [g.a*, ]= 0 


where, |®) is the ground state and: 


AN, 


[n,], - (16) 
—&§ 


[n, ],!= [n, |, [n, =A ee [E (17) 
The eigenstates I") of the oscillator number 
operators is as follow: 
N,|n) = n,|n) (18) 
From Eq. (14), it can be seen that the action of 
operators 4; and “* on the eigenstates |”) gives: 
a; |n) = fln.+1],|7, +1,n_,) 
a,|n)= lr], —1) (19) 
a’,|n)= 8" Vest, |mor +1) 
a.g|n)=8" lnc], |eon-4 -1) 


The energy spectrum of &-deformed crystal lattice 


vibration for a generic atomic string is given by: 
H\n) =E, n) (20) 


where, 
1 


E, =) AO fn, +1), +l], } (21) 


2 


If the Hamitonian (13) of the crystal lattice vibration 
for a generic atomic string in the &-deformed quantum 
theory can be expressed in terms of creation 
annihilation operators 7; , 4%, which are defined to 
satisfy the commutation relations (14), the crystal 
lattice vibration are then considered as particle system 


and its energy spectrum depended on deformed 
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parameters according to Eq. (21). 


4. The $-Deformed Green Function 


The g-deformed Green function is defined as the 
statistical distribution of 4:4: : 


+ 1 -BoN, ~+ 
(aja, )= Str (e "apa, ) (22) 
where, Z is the partition function: 


Z=Tr (err )= y (n lereMe 


n=0 


n) (23) 


which determines the thermodynamic properties of the 
system, £ = I/kT, w being one particle-oscillator 
energy. The trace must be taken over a complete set of 
states. 

The caculations based on the Eq. (15) give the 
following results: 


1 
pe (24) 


5. Conclusions 


The present construction a statistical distribution 
for &-deformed crystal lattice vibration for generic 


atomic string is an example of use deformed 
oscillators formalism to study of condensed matter 
physics, this is a fundament for study of crystal lattice 
vibration, which would give a base for a new 
quantization procedure. Further the next step is to 
extend the present investigations to the case of other 
physical systems. 
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